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1.0 INTRODUCTION 

BOPACE is the acronym for the Boeing plastic Analysis £apability for 
Engines. BOPACE was developed by Boeing/Huntsville to meet the evident 
need for an advanced thermal -el as tic-plastic-creep structural analyzer. 
Although BOPACE development has been strongly influenced by the require- 
ments for structural analysis of engines, in particular the space shuttle 
main engine, its capabilities have been kept quite general and it is 
applicable to many types of nonlinear structures. 

The philosophy for program development was based on the following 
requ i remen ts . 

1) Analysis of very high temperature and large plastic-creep effects. 

2) Treatment of cyclic thermal and mechanical loads. 

3) Improved material constitutive theory which closely follows actual 
behavior under variable temperature conditions. 

4) A stable numerical solution approach which avoids cumulative errors. 

5) Capability for handling up to 1000 degrees of freedom with moderate 
computation cost. 

Although the finite-element method was first applied to plasticity 
in the early 1960's, and several good programs for nonlinear analysis 
have since been developed, numerous improvements were indicated in order 
to satisfy the above requirements. For example, some other available 
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1.0 (Continued) 

programs assume linear plastic hardening, accumulate errors by failing 
to satisfy equilibrium at each step, or do not completely account for 
the effects of variable temperature on the elastic and plastic relations. 
The stated requirements have been effectively met by the current BOPACE 
program version. In addition, the research and development effort has 
led to an improved hardening theory for cyclic plasticity, a method for 
representing general cases of load reversal, and advanced techniques for 
improving the accuracy and controlling convergence of highly nonlinear 
solutions. 

Two versions of the current BOPACE program are available. The first 
is a 300-D0F version developed for fast analysis of small size problems 
within moderate core-storage limitations. The second is the basic 
lOOO-DOF version. In addition, a low-core modification of the lOOO-DOF 
version has been accomplished through the use of overlays and dynamic 
storage of arrays. BOPACE is written in FORTRAN IV and has been 
extensively run on both the IBM 360 and UN I VAC 1108 computer systems. 
Documentation consists of three volumes: Theoretical Manual, User 

Manual (including example problems), and Programmer Manual. 

The BOPACE development and programming effort has been performed at 
Boeing/Huntsville by Dr. R. G. Vos, with suggestions and review by 
W. H. Armstrong. A. H. Spring assisted with many analyses and program 
checkout. 0. L. Ballinger of Boeing Computer Services modified and 
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1.0 (Continued) 

programned the Gauss wavefront solution method. Recognition is also due 
to N. L. Schlemmer, L. Salter and R. Hurford at the NASA Marshall Space 
FLight Center, and L. Johnston of Brown Engineering Co., for their 
suggestions and support of the program development. 
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2.0 SUMMARY OF BOPACE IMPUT DATA 

A pictorial of the BOPACE input deck is shown in Figure 2.0-1. The 
input data consists of the following three general types: 

Type C: Data on the usual card file. These are data 

which are needed for each start or restart. 

Data on File I. These are basic structural 
data for a given problem, such as material 
properties and mesh data. They are the same 
for all load increments and are needed only 
when starting. 

Data on File II. These are incremental thermal - 
load and z-load data which are needed for each 
start or restart. 

The data included on each file are described below. Formats are 
consistent with FORTRAN IV conventions, 

C-1 . Start-restart code and data file numbers: 

a. “START" if new problem, or "RESTART" if restarting. 

b. If starting give unit number for file I. 

c. Unit number for file II. 

d. Unit number for output file (e.g. printer). 


Type I: 


Type II: 
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2.0 (Continued) 

e. If restarting give load increment number from the end of which 
a restart is to be made. 

f. If restarting give input restart- tape unit number. 

g. If data is to be saved for future restart give output restart- 
tape unit number. 

Format (A4,6X,6I5) 

C-2. Problem I.D. title. 

Format (20A4). 

C-3. Program control constants (any constant left blank is assigned 
a default value): 

a. Code for system matrix decomposition and solution. This 
code controls only the method of iteration and convergence, 
and does not affect final computed results. 

Code 1 = use only elastic matrix with no updating. 

2 = update elastic matrix. 

3 = update plastic matrix. 

4 - update total Jacobian (elastic+pl astic) matrix. 

5 = update both elastic and total Jacobian matrices. 

The default code is 5. 
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1-7 ELEHENT DEFINITIONS 


1-6 NODE DEFINITIONS 



C-5 FINAL BLANK CARO 

II-l THERMAL AND Z LOADS FOR 
EACH INCREMENT 


C-4 MECHANICAL LOAD DATA 


1-9 MECHANICAL LOAD 
REFERENCE VECTORS 


1-8 FORCE -DISPLACEMENT- 

CONSTRAINT SPECIFICATIONS 


1-5 SPECIAL COORDINATES 


1-4 CREEP DATA 


1-3 PLASTICITY DATA 


1-2 WTERIAL PROPERTIES 


I-l STRUCTURE COOES 
AND CONSTMTTS 

C-3 CONTROL CONSTANTS 


C-2 PROeiEM I.D. 

C-1 START-RESTART 

CODE AND FILE NMCERS 


~h 





NOTES: 

C IS CARD FILE. 

I. II ARE CARD OR TAPE FILES. 

DATA ON FILE I IS USED ONLY 
WEN STARTING. DATA ON FILE II 
IS USED WEN STARTING OR 
RESTARTING. 


FIGURE 2.0-1: BOPACE INPUT DECK SETUP 
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2.0 (Continued) 

b. Maximum number of stiffness updates per load increment (in 
order to achieve convergence to within allowable error norm). 
Defaul t = 1 . 

c. Maximum number of residual -force iterations per stiffness 
update. Default = 10. 

d. Maximum number of initial iterations using elastic matrix 
(to account for possible unloading). Default = 2. 

e. Maximum allowable magnitude for elastic-plastic sum code 
YCODEl. Default = 2. 

f. Maximum number of cuts to be performed (giving new solution 
with a fraction of previously used displacement corrections) 
if error norm is not decreasing. Default = 0. 

g. Cutting fraction (displacement correction = previous correction 
times cutting fraction). Default = 0.5. 

. h. Maximum allowable error norm. Default = 0.001. 

i. Fraction from end of increment to evaluate stress vs. plastic- 
strain slope in forming plastic stiffness. Default = 0.1. 

Format (6I5,3F10.0) 

I-l. Plane-stress code (0) or plane-strain code (1), number of materials, 

default thickness, fabrication temperature. 

Format (2I5,2F1Q.0) 
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2.0 (Continued) 

1-2. Material property data for each consecutive material. 

a. Material number 
Format (110) 

b. Three consecutive temperature-dependent property curves 
(thermal strain, elastic modulus, Poisson's ratio). For 
each curve point give temperature and value, with points 
in order of increasing temperature. User has option of 
from 1 to 4 points per card. 

Format (8F10.0) 

Blank card after last point of each curve. 

1-3. Plasticity data for each consecutive material. 

a. Material number, plasticity type, kinematic code. 

Type 1 = strain hardening 
2 = work hardening 

Code 0 = kinematic hardening is function of one parameter 
1 = kinematic hardening is function of two parameters 

Format (3110) 

b. Temperature-parameter-hardening data, in order of low to high 
temperatures. For each temperature: 
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2.0 (Continued) 

Material number, temperature. 

Format (IIO.FIO.O) 

Three consecutive hardening curves (cumulative hardening 
parameter vs. isotropic hardening, i.e. yield surface size; 
kinematic parameter vs. kinematic hardening shape; cumulative 
parameter vs. kinematic hardening factor). Bauschinger 
hardening is computed as kinematic hardening shape times 
hardening factor. If kinematic code = 0, curve of kinematic 
factors is not given and all factors are taken as 1.0. For 
each point give parameter and hardening value. User has 
option of from 1 to 4 points per card. First point on each 
curve must be yield point (parameter = 0.0). 

Format (8F10.0) 

Insert blank card after each input curve. 

Blank card after all temperatures for a given material. 

1-4. ' Creep data for each consecutive material. 

a. Material number, creep type. 

Type 1 = age hardening 

2 = strain hardening 

3 - work hardening 
Format (2110) 
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2.0 (Continued) 

b. Reference creep curve. For each point, the time and creep 
strain, in order of increasing time. User has the option 
of from 1 to 4 points per card. 

Format (8F10.0) 

Blank card after last point of curve. 

c. Table of creep factors, in order of low to high temperatures. 
Creep is computed as factor (function of temperature and 
stress) times reference creep curve. 

For each temperature: 

Material number, temperature 
Format (IIO.FIO.O) 

For each point given at this temperature, the stress 
and creep factor, in order of increasing stress. User 
has option of from 1 to 4 points per card. 

Format (8F10.0) 

Blank card after all points for a given temperature. 

Blank card after all temperatures for a given material. 

1-5. For each special Cartesian coordinate system; the identification 
number (integer = 2) and counter-clockwise angle (degrees) from 
basic system X-axis to special system x~axis. 

Format (I10,F10.0) 
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2.0 (Continued) 

Blank card after last coordinate system. 

1-6. For each node: Node I.D. number, identification number of 

coordinate system to define location, X and Y (or R and 6), 
identification number of coordinate system to define displacements. 
(Coordinate identification number 0 implies the basic Cartesian 
system, 1 implies the basic cylindrical system). Order of nodes 
in data deck is internal order used to form system matrix. 

Format (2I5,2F10.0, 15) 

Blank card after last node. 

1-7. For each element: element I.D. number, material number, thickness, 

three node I.D. numbers (counter-clockwise order). If thickness 
is left blank, default value from I-l is used. 

Format (2I5,F10.0,3I5) . 

Blank card after last element. 

1-8. . For each degree of freedom with a specified force, displacement 
or constraint: give node I.D. number, component number (1 or 2), 
and code. The code to be given is: 

1. For specified force, the node I.D. number 

2. For specified displacement, the negative of the 
node I.D. number. 
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2,0 (Continued) 

3. For dependent constrained DOF, the node I.D. 
number of the independent DOF in constraint. 

The default code is specified force. User has option of from one to 
four DOF per card. 

Format (4(315, 5X)) 

Blank card after last force-displacement-constraint DOF. 

1-9. Mechanical load reference vectors 

Number of vectors (for current program version must be 2) 

Format (110) 

For each non-zero component of load vector: node I.D. number, 

component number (1 = X or R, 2 = Y or 6), value. User has 
option of from 1 to 4 values per card. 

Format (4(215, FIO.O)) 

Blank card after last value of each vector. 

C-4. Incremental mechanical load data. 

Number of load increments 
Format (110) 

For each load increment: maximum iterations per stiffness update 

(if left blank, value from C-3 is used), the cumulative factors 
to be applied to load reference vectors (for current version of 
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2.0 (Continued) 

program two factors must be given), creep time increment (if left 
blank, no creep calculations are made). 

Format (I10,3F10.0) 

II-l. Incremental thermal and z-di recti on load data. 

a. Increment I.D. title 
Format (20A4) 

b. Element thermal loads. For each specified component of load: 
element I.D. number, temperature at end of increment. Non- 
specified temperatures are taken to be the fabrication 
temperature for the first increment; for later increments they 
are taken to be the temperature of the preceding increment. 
User has option of from one to four values per card. 

Format (4(110, FIO. 0) ) 

Blank card after last specified value of thermal load for each 
load increment. 

c. Element z-stress (for a plane-stress problem) or z-strain (for 

a plane-strain problem). Plane-stress or plane-strain condi- 
tion is defined by code in data item I-l. For each specified 
element z-load: element I.D. number, stress or strain at the 

end of increment. Non-specified z-loads are taken to be zero 
for the first increnent; for later incronents they are taken 
to be the z-load of the preceding increment. User has option 
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2.0 (Continued) 

of from one to four values per card. 

Format (4(110, FlO.O)) 

Blank card after last specified value of z-load for each load 
increment. 


C-5. Blank card after last problem. 
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3.0 SUMMARY OF OUTPUT 

A discussion of BOPACE output is conveniently divided into two parts. The 
first covers output which is primarily an echo check of the input data, 
and the second covers output results for each load increment. 

3.1 ECHO CHECK OF INPUT DATA 

Initial Output - The first page of BOPACE output for a problem is es- 
sentially an echo check of input items C-1, C-2, C-3 and I-l. An indica- 
tion is given as to whether the problem is being started or restarted. If 
it is restarted then the previous increment number is given, from the end 
of which the restart is progressing. Next the problem I.D. title is 
printed, followed by the various control constants which are determined 
from their default values or from input values C-3. Finally the input 
data from I-l are printed. 

Material Properties - For each material the three curves (thermal strain, 
elastic modulus and Poisson's ratio) input in data item 1-2 are printed. 

Plasticity Data - For each material the plasticity data input in data item 
1-3 are printed. These include the plasticity hardening type (1 = strain 
hardening, 2 = work hardening) and the kinematic hardening code (0, 1 for 
one, two parameter hardening, respectively). Following are groups of two 
or three curves given at each temperature (surface size a - a vs. cumula- 
tive hardening parameter <, kinematic hardening curve shape vs. kinematic 
shape parameter k , and kinematic curve magnitude vs. k if given). The 
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3.1 (Continued) 

hardening curves are input and output by order of associated low to high 
temperature. Abscissas of the curves in the first (lowest temperature) 
group are used as the basis for tabulating all curves* i.e., curve points 
of higher temperatures are interpolated to the low-temperature abscissas. 
These interpolated hardening curves are printed following printing of 
the input hardening curves. 

Creep Data - For each material the creep data input in data item 1-4 are 
printed. These include the creep hardening type (1, 2, 3 for age, strain, 
work hardening, respectively), and the reference creep curve shape of 
creep strain vs. time. Following are the creep factors for each combina- 
tion of temperature and stress, grouped by low to high temperature. Stress 
values in the first (lowest temperature) group are used as the basis for 
tabulating all groups, i.e., creep factors at higher temperatures are 
interpolated to the low-temperature stress values. The interpolated harden- 
ing factors are printed following printing of the input hardening factors. 

Special Coordinate Systems - These are the user-defined direction (special 
Cartesian) systems of input data item 1-5. Quantities printed are the 
system I.D. number, and counter-clockwise angle (in degrees) from the 
basic X axis to the special -system x axis. 

Node Definitions - The information given in input item 1-6 is printed. 

Values are the node number, node I.D. , location coordinate system number 
(0 = basic Cartesian, 1 = basic cylindrical), X or R coordinate, Y or 0 
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3.1 (Continued) 

(degrees) coordinate, direction coordinate system number (0 = basic 
Cartesian, 1 = basic cylindrical, >1 = I.O. of special system). 

Element Definitions - The information given in input item 1-7 is printed. 
Values are the element number, element I.D. , material number, element 
thickness, the three element node I.D. numbers in counter-clockwise order, 
and the computed element area. 

Force-Displacement-Constraint Prescriptions - These are the codes given 
in input data item 1-8. Quantities printed are the node I.D., the com- 
ponent nuni)er, and code for each degree of freedom with a user-specified 
force, displacement or constraint. 

Mechanical Load Reference Vectors - For each input component of the two 
load vectors from input item 1-9, the node I.D. , component number, and load 
are printed. 

Incremental Mechanical Load Data - Quantities related to input data item 
C-4 are printed. First is printed the number of load increments to be run. 
Then for each increment is given the increment number, input or default 
value for maximum number of iterations per increment, factors to be applied 
to the two load reference vectors, and the creep time increment. 
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3.2 RESULTS FOR EACH LOAD INCREMENT 

Thermal and Z Loads - The cumulative values of thermal load and 1 - 
direction mechanical load for each element, given in input item 11*1, are 
printed. A heading gives the increment number, and is followed by the 
input I.D. title for the increment. Then the element I.D. , cumulative 
temperature and cumulative Z load are printed, by groups of ten elements. 

Iterative Error Values - An error norm computed at the end of each 
iteration is printed. The error norm is obtained by a ratio of unbalanced 
(residual) forces to "total” forces. The total forces are computed by 
summing a stress quantity times the thickness of each element, multiplied 
by the square root of the element area. This gives a meaningful error 
norm even if the applied loads do not include any forces. 

Increment Heading - The load increment number is printed, along with the 
increment I.D. title given in input item II-l. Following this are the 
creep time increment, the number of elastic and plastic elements at the 
end of the increment, the number of elements which have changed elastic 
to plastic and plastic to elastic during the increment, the maximum allowable 
number of Jacobian updates and the number performed during this increment, 
the maximum allowable number of iterations per update and the number per- 
formed since the last update, and the maximum allowable error norm and 
the error norm actually obtained. 

Cumulative Forces and Displacements - These are the cumulative internal 
forces (corresponding to computed element stresses) and displacements for 
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3.2 (Continued) 

each node. The node number and node I.D. are printed, followed by the 
U and V components of force and displacement. U and V are in the directions 
defined by the node direction coordinate system (X-Y, R-9, or special 
coordinate system). 

Thermal and Elastic Strains - These are incremental and cumulative values 
of the thermal and elastic strains for each element. The element number 
and element I.D. are printed, followed by the thermal strains and components 
of the elastic strains. All strain components are referred to the element 
base coordinate system (an x-y right-hand Cartesian system with origin at 
node 1, and x axis along the nodal line 1-2). 

Plasticity Results - For each element the number and I.D. are printed, 
followed by the incremental and cumulative plastic work and components of 
plastic strain. All strains are again given in the element base 
coordinate system. 

Creep Results - These are printed only if the creep time increment 
(input in data item C-4) is greater than zero. Printed creep results 
correspond to those for plasticity. 

Stress Results - For each element the number and I.D. are printed, 
followed by the values of cumulative effective stress center a 
and effective stress o, and the components of stress center and stress 
referred to the element base coordinate system. 
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3.2 (Continued) 

Summarized Element Quantities - For each element the number and I.D. 
are printed, followed by the elastic-plastic "code" and "sum code." 

The code is 0, -1 or +1, respectively, according to whether the element 
condition has remained unchanged, gone from plastic to elastic (unloaded), 
or gone from elastic to plastic (yielded) during the increment. The 
sum code gives the value for the program variable YCODEl, and is + or 
respectively, according to whether the element condition is plastic or 
elastic at the end of the increment. Its magnitude is an indication of 
the iterative tendency for the element to remain in that condition. 

Mext are given the total temperature and the yield surface size (a - a) 
at the end of the increment. Finally are printed, for both plastic and 
creep strains, three values of effective strain (incremental aF, sum of 
incremental lAe, and cumulative e). 
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4.0 SIZE LIMITATIONS 

4.1 GENERAL SIZE LIMITATIONS 


The following variables are used to specify maximum size limitations in 
BOPACE. The values set for these variables in the 300-D0F and lOOO-DOF 
program versions are given in Table 4.1-1. 


NSTOR = 
NMAXl 
NMAX2 = 
NMAX3 = 
NMAX4 = 
NMAX5 = 
NMAX6 = 
NMAX7 = 

NMAX8A = 
NMAX8B = 

NMAX8C = 

NMAX9 = 
NMAXl 0 = 
NMAXl 1 = 

NMAXl 2 = 
NMAXl 3 = 


core words allocated to Gauss wavefront storage 

maximum number of materials 

maximum number of nodes 

maximum number of elements 

maximum node I.D. number 

maximum element I.D. number 

maximum number of points in a material property curve 

maximum number of temperature plasticity curves per 

material 

maximum number of points per isotropic hardening curve 

maximum number of points per kinematic hardening shape 

curve 

maximum number of points per kinematic hardening factor 
curve 

maximum number of points in a creep reference curve 
maximum number of creep-factor temperatures per material 
maximum number of creep-factor stresses per temperature 
maximum number of special coordinate systems 
maximum (required) number of mechanical load reference 
vectors 
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4.1 (Continued) 

NMAX14 = maximum number of load increments per run 


TABLE 4.1-1: 

MAXIMUM SIZE 

LIMITATIONS 

Maximum 

300 DOF 

1000 DOF 

NSTOR 

2000 

5000 

NMAXl 

5 

5 

NMAX2 

150 

500 

NMAX3 

200 

800 

NMAX4 

500 

2000 

NMAX5 

1000 

3000 

NMAX6 

20 

20 

NMAX7 

6 

6 

NMAX8A 

30 

30 

NMAX8B 

20 

20 

NMAX8C 

30 

30 

NMAX9 

10 

10 

NMAXIO 

6 

6 

NMAXl 1 

10 

10 

NMAXl 2 

50 

50 

NMAXl 3 

2 

2 

NMAXl 4 

60 

60 
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4.2 LINEAR EQUATION SOLVER LIMITATIONS 

The linear equations solution routines have one user controlled limi- 
tation. This is the maximum bandwidth of active nodes during the decom- 
position of the stiffness matrix. The bandwidth is defined as the number 
of nodes following the node being processed which have non -zero terms 
associated with the node being processed. Melosh and Bamford [1] discuss 
this in some detail as the wavefront analysis concept. Whetstone [2] 
also discusses this concept and gives rules and procedures which can be 
used to keep the bandwidth as small as possible by proper numbering of the 
nodes . 

BOPACE uses the maximum active bandwidth (wavefront) to determine core 
storage requirements during decomposition. 

Let 

K = the amount of storage available (K = 2000 for the 300 DOF 
version, K = 5000 for the 1000 DOF version) 

N = the number of DOF per node (a constant for each node in 

the analysis) 

B = the maximum active bandwidth 

T = B+B*N^+4 

M = B * (B+D/2 
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4.2 (Continued) 

Then the following two equations must be satisfied: 
T < 420 

K > B + M + 2*N^*M + 6*N^ + T 


(4.2-1) 

(4.2-2) 
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5.0 DETAILED INPUT AND STORAGE OF DATA 

5J DETAILED INPUT DESCRIPTION 

This section provides a definition of BOPACE input variables, and a de- 
tailed description of the input data and their use within the program. 

Data are discussed in the order in which they are read, with data item 
numbers corresponding to those given in Section 2.0. 

Input Files - In addition to the usual card file (unit number 5), BOPACE 
uses files I and II for which the user defines unit numbers. Because the 
data on these files may become rather lengthy for a large-size problem, 
the user may wish to define them as tape files which are automatically 
generated by special subroutines. (See for example Appendix A, for a 
description of the INPUTS interpol ation-data-generator program.) File I 
data define the basic problem, mesh and material properties. These data 
cannot be redefined during solution of the problem, so that they are in- 
dependent of any particular load increment and are used only when starting 
a new problsn. 

C-1. Start-Restart Code and Data File Numbers - These data are read by 
subroutine READRS from the usual card file (unit number 5): 


a. 

START 

= 

start-restart code 

b. 

UINl 

= 

unit number for input file I 

c. 

UIN2 


unit number for input file II 

d. 

UOUT 

= 

unit number for output file (e.g., printer) 
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5.1 (Continued) 

e. INCR = previous increment number from the end of which 

a restart is to be made 

f. UINRS = unit number for input restart tape 

g. UOUTRS = unit number for output restart tape 

A "RESTART" code specified in item (a) means that an input restart tape is 
being provided (item f), which contains the initial Jacobian and its 
decomposition, the basic problem, mesh and material data, and previous 
incremental results including those for increment INCR (e) from which a 
restart is to be made. If INCR = 0, the problem is started from the 
initial conditions, but the reading of the file I data and the formation 
and decomposition of the initial Jacobian are avoided. If neither "START" 
nor "RESTART" is specified, it is assumed that the blank card of input 
item C-5 (see Section 2.0) has been read, and the program exits in the 
normal mode with a STOP 9999 code. 

The file unit number (b) is required when starting. Unit numbers (c) and 
(d) are required in all cases. If the current run is a restart, the unit 
number (f) is required. If results are to be saved for a future restart, 
the unit number (g) is required. 

C-2. Problem I.D. Title - The title is read by subroutine READO: 

I DENT = problem I.D. title 

The title consists of any 80 characters used to describe the current run. 
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C-3. Program Control Constants - These constants are read by subroutine 
READO: 


a. 

SCODE 

= 

code for system matrix (Jacobian) updating 

b. 

MAX UP 

= 

maxi mum 

number of updates per increment 

c. 

MAX IT 

= 

maximum 

number of iterations per update 

d. 

MAX IE 

= 

maximum 

number of initial elastic iterations 

a 

MAXYC 

= 

maximum 

allowable magnitude for YCODEl 

f. 

MAXCUT 

= 

maximum 

number of cuts 

g- 

CUT 

= 

cutting 

fraction 

h. 

ERRMAX 

= 

maximum 

allowable error norm 

1 . 

AFACT 

= 

Increment fraction for evaluating hardening slope 


The value of SCODE (item a) is the code for updating the Jacobian and its 
component matrices. It allows various options for the basic solution 
approach, depending on the amount of nonlinearity expected in the problem. 

It affects only the rate of convergence and the effectiveness of the 
iterative process. It does not alter the manner in which unbalanced 
forces are iteratively computed, nor does it affect the final computed 
results so long as convergence is achieved. A discussion of the effects 
of each option for various types of problems is given in Section 4.4 of the 
BOPACE Theoretical Manual. Item (b) gives the maximum number of Jacobian 
updates per load increment. An update is performed if the maximum number 
of iterations specified in (c) has been reached without achieving convergence. 
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5.1 (Continued) 

Generally a maximum of one update per increment is sufficient, unless 
very large changes in the elastic properties, slope of the stress-strain 
curve, or direction of loading, have occurred. 

Item (d) gives the maximum number of iterations to be performed at the 
start of the increment, using the current elastic matrix. These initial 
iterations are performed in order to account for the possible occurrence 
of large-scale unloading within the structure. 

MAXYC (item e) gives the maximum absolute value for the elastic-plastic 
sum code, YCODEl. The current value of YCODEl for each element deter- 
mines whether it is to be treated as elastic or plastic in the calcula- 
tion of unbalanced forces and in the formation of the Jacobian stiffness 
matrix (0 or - denotes elastic condition, + denotes plastic condition). 
YCODEl is changed by -1 or +1 after each iteration, depending on 
whether the condition was determined to be elastic or plastic, respectively, 
but its absolute value is not allowed to exceed MAXYC. If the value of 
YCODEl reaches 0, the element condition is considered to have changed 
and YCODEl is set to an initial value of ±2 (±2 is used instead of ±1 
because it prevents oscillation between elastic and plastic conditions 
in certain cases). The purpose of YCODEl is essentially to denote the 
element condition and to provide a stabilizing influence on the iterative 
process. Higher specified values of MAXYC could increase the stability 
but might require larger convergence times. 
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Items (f) and (g) (MAXCUT and CUT) relate to a common method for 
increasing the stability of iterative processes. If a computed displacement 
correction in the BOPACE iteration does not result in a decreased error 
norm, then the error norm is computed for another configuration based 
on some fraction of the displacement correction. If the error norm still 
has not decreased, the fraction is squared and again applied to the 
displacement correction, etc. MAXCUT gives the maximum number of such 
"cuts” to be performed, and CUT gives the fraction to be used. 

Item (h) gives the maximum allowable error norm. The actual error norm 
is computed essentially as a ratio of unbalanced forces to applied forces. 

AFACT (item i) relates to evaluation of the plastic hardening slope for 
each element, used in forming the Jacobian matrix. The best location for 
evaluating this slope for the Jacobian update is a point at the end of the 
current increment. The actual slope is determined in BOPACE using two 
interpolated values from the hardening tables. The second corresponds 
to the estimated value of the hardening parameter at the end of the 
increment, while the first is at a small distance prior to the second and 
is determined by subtracting the fraction AFACT times the incremental 
change in the hardening parameter. 

I-l. Basic Constants and Default Values - These items are read by 
subroutine READl : 
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5.1 (Continued) 


PCODE = 

plane-stress or plane-strain code 

NMAT 

number of materials 

THICK = 

default thickness 

TEMPO = 

fabrication temperature 


The code for plane-stress or plane-strain holds for all elements and 
defines the type of problem being run. The default thickness defines 
the thickness of all elements, unless it is overidden by the element 
values given in 1-7. The fabrication temperature defines the initial 
temperature of the structure, at which thermal strains are considered 
to be zero. It is also the default temperature for all elements during 
the first load increment, unless it is overridden by the element values 
given in Il-lb. 

1-2. Material Property Data - These data are read by subroutine READTM: 


NTHERM(I) 

number of points in thermal-strain curve for 
material I 

THERMX(0,I) = 

Jth temperature value of thermal -strain curve for 
material I 

THERMY(J,I) = 

Jth thermal strain value for material I 

NEMOD(I) 

number of points in elastic-modulus curve for 
material I 
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EM0DX(J,I) = Jth temperature value of elastic-modulus curve 

for material I 

EM0DY(J,I) = Jth elastic-modulus value for material I 

NPRAT(I) = number of points in Poisson's ratio curve for 

material I 

PRATX(J,I) = Jth temperature value of Poisson's ratio curve 

for material I 

PRATY(J,I) = Jth Poisson's ratio value for material I 

The property data are given consecutively for each material. The zero- 
strain point of the thermal -strain curve is arbitrary, because the curve 
is used only to obtain the increment of thermal strain corresponding to 
a given temperature increment. All interpolation performed later from 
the thermal -strain, elastic-modulus and Poisson's ratio curves is 
accomplished by the linear interpolation routine YVAL. 

1-3. Plasticity Data - The plasticity data are given consecutively 
for each material. The data are read by subroutine READTP: 

PTYPE(I) = plastic hardening type for material I 

KTYPE(I) = kinematic hardening code for material I 

NTABY(I) = number of temperature values for hardening 

tables for material I 
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TABY(J,I) = Jth temperature value for hardening tables for 

material I 

NITABX(I) = number of k values for isotropic hardening table 

for material I 

ITABX(J,I) = Jth K value for isotropic hardening table for 
material I 

ISTAB(K,J,I) = value of isotropic hardening table for Kth tem- 
perature value and Jth k value for material I 

NKTABX(I) = number of k*^ values for kinematic hardening 

shape table for material I 

KTABX(J>I) = Jth value for kinematic hardening shape table 
for material I 

KSTAB(K,J,I) = value of kinematic hardening shape table for Kth 
temperature value and Jth value for material I 

NFTABX(I) = number of k values for kinematic hardening factor 

table for material I 

FTABX(J,I) = Jth < value for kinematic hardening factor table 
for material I 

FSTAB(K,J,I) = value of kinematic hardening factor table for Kth 
temperature value and Jth k value for material I. 
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PTYPE defines for each material whether it is a strain-hardening or work- 
hardening material. KTYPE defines whether the kinematic hardening for 
each material is to be taken as a function of only the parameter k^, or 
of both and k. Two-parameter kinematic hardening accounts for change 
in magnitude of the Bauschinger effect with continued cycling, i.e., 
with increase in value of k. 

Hardening is defined by values stored in the three hardening tables 
ISTAB, KSTAB and FSTAB. ISTAB gives the size of the yield surface 
(a - a). KSTAB gives the shape of the kinematic hardening (a) curve, 
and FSTAB gives the factor to be applied to this shape in order to 
define the variable Bauschinger effect. The values in these 2-dimensional 
tables (for each material) represent a series of hardening curves at 
various temperature levels. TABY defines the temperature series which 
is used for all three hardening tables. However, each table has its own 
set of values for hardening parameter, defined by the respective vectors 
ITABX, KTABX and FTABX for each material. When, for example, the k- 
parameter ITABX and the temperature value TABY are given for a particular 
material, the corresponding size of the yield surface is determined by 
double interpolation from the table ISTAB. 

The hardening curves are input in order of associated low to high 
temperature. The curve at each temperature may have its own series of 
K or K parameters. However, a single series of these parameters (defined 
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by ITABX, KTABX or FTABX) is needed for each hardening table. This 
series is taken as that given for the low- temperature curve of the table, 
and curves for other temperatures are interpolated to these basic values. 
Therefore the low- temperature curve input for each table should cover 
the range of parameter values used in all curves input for that table, 
in order to avoid losing any points in the initial interpolation and forming 
of the table. Later interpolation from the plastic hardening tables is 
accomplished by the double linear interpolation routine ZVAL. 

For KTYPE = 0, kinematic hardening is determined using only the KSTAB 
table, and the FSTAB table and associated data are not read. For KTYPE = 1, 
kinematic hardening is determined as a product of values from the shape 
table KSTAB and the factor table FSTAB. A more detailed description of 
the BOPACE plastic hardening method is given in Section 2.3 of the 
Theoretical Manual. 


1-4. Creep Data - The creep data are given consecutively for each material. 
The data are read by subroutine READTC: 


CTYPE{I) 
NCREEP(I) = 


CREEPX(J,I) = 
CREEPY(J,I) = 


creep hardening type for material I 
number of points in reference creep curve for 
material I 

Jth time value of reference creep curve for material I 
Jth creep-strain value of reference creep curve for 
material I 
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NCTABX(I) 

number of stress values for creep factor table 
for material I 

NCTABY(I) 

number of temperature values for creep factor 
table for material I 

CTABX(J,I) 

Jth stress value for creep factor table for 
material I 

CTABY(J,I) 

Jth temperature value for creep factor table 
for material I 

CTAB(K,J,I) = 

value of creep factor table for Kth temperature 
value and Jth stress value for material I 


CTYPE defines for each material whether it is an age-hardening, strain- 
hardening or work-hardening material with respect to creep. 

Shape of the creep curve is defined by the vectors CREEPX and CREEPY 
for each material. Creep factors are stored in the 2-dimensional table 
CTAB. CTABY defines the temperatures associated with these factors, 
and CTABX defines the stress levels. Later interpolation from the CTAB 
factor table is accomplished by the double linear interpolation routine 
ZVAL. Interpolation from the shape vector CREEPY is accomplished by 
the linear incremental interpolation routine DYVAL, which gives the 
increment of strain on the reference creep curve corresponding to given 
initial and final values of time. The actual increment of effective 
creep strain is determined as a product of the CTAB factor and the CREEPY 
shape increment. 
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Each input creep hardening factor has an associated temperature and 
stress level. The factors are input in temperature groups, in order 
of low to high temperature, and within each group are ordered from low 
to high stress level. Although each temperature group may have its own 
input stress-level series, the stress-level series defined by CTABX is 
used for all temperatures during later interpolation. CTABX is taken 
as the stress series given for the low- temperature group of factors, 
and other groups are interpolated to these basic values. Thus the input 
low- temperature group should cover the range of stress values input 
for all groups, in order to avoid losing any points in the initial 
interpolation and forming of the factor table. 

1-5. Special Coordinate Systems - The special coordinate systems are 
read by subroutine READC: 

COORDA(I) = angle for special coordinate system I 

These systems provide special nodal direction coordinates which can be 
used to specify components of force or displacement in a particular 
direction. 

1-6. Nodal Data - These data are read and computed by subroutine READM: 

NOD = number of nodes in problem 

COORD(J,I) = Jth (X or Y) coordinate for node I 

GC0S(J,I) = Jth direction Cosine for node I 
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NODE(I) = node I.D. for node I 

NODI(I) = node number for node I.D. I 

The value of NOD is determined by the program from counting the input 
nodal data. When the two coordinates for a node are read, they are con- 
verted to basic Cartesian (X-Y) coordinates and stored in the COORD 
array. The GCOS direction Cosines are the Cosine and Sine, respectively, 
of the angle from the basic X axis to the nodal direction coordinate 
system X axis. The GCOS array is computed and stored in double 
precision, so as to avoid loss of accuracy. This precision is required 
in problems such as those involving a small cylindrical segment of an 
engine wall, where quantities involving the Sine and Cosine of very 
small angles are important. 

The node I.D. is supplied by the user, and actual internal node numbers 
are assigned consecutively in the order in which nodes are read in the 
input data deck. Since the internal numbers are used in all operations 
including formation and solution of the linear equations, the user 
should attempt to place the nodes in the data deck in an efficient order, 
i.e., an order which avoids as much as possible the connection of nodes 
having greatly different internal node numbers. 

1-7. Element Data - These data are read and computed by subroutine READM: 


NEL 

IMAT(I) 


number of elements in problem 

material number for element I 
5-13 



05 ^ 17266^2 


5.1 (Continued) 


T(I) 

ELN0(J,I) = 
NELE(I) = 
NELI(I) = 


thickness of element I 
Jth node number for element I 
element I.D. for element I 
element number for element I.D. I 


The value of NEL is determined by the program from counting the input 

elanent data. The ELNO array stores the three internal node numbers for 

/ 

each element in counter-clockwise order. 

The element I.D. is supplied by the user* and actual internal element 
numbers are assigned consecutively in the same manner as for nodes. In 
order to increase the efficiency of the global matrix formations, it is 
best to have the order of elements in the data deck follow approximately 
the order of the nodes to which they are connected. 


1-8. Force-Displacement-Constraint Specifications - These data are 
read by subroutine READ2: 

KFD(I) = force-displacement-constraint code for degree of 
freedom I 

The degrees of freedom are ordered by node number, e.g., KFD(5) is the 
code for the first degree of freedom at node 3 (2 DOF per node). Each 
DOF has either a specified force or displacement. For the usual type 
of structure, most will be specified forces. For a constrained DOF, 
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i.e., a dependent DOF whose displacement is set equal to that of some 
independent DOF, the program makes its code agree with that of the 
independent DOF. Thus if the independent DOF has a specified force, the 
dependent DOF will have a specified (perhaps zero value of) force; 
if the independent DOF has a specified displacement, the dependent DOF 
by definition has the same displacement. 

A DOF may not be constrained with a dependent (already constrained) DOF. 
Also, the present BOPACE equation-solver routines do not allow constraints 
between two nodes on the same element. 

1-9. Mechanical Load Reference Vectors - These data are read by sub- 
routine READS: 

PREF(J,I) = load component of Jth (1st or 2nd) load reference 
vector for DOF I 

The load components are given in the directions defined by the node 
direction coordinate systems. The current version of BOPACE employs two 
load reference vectors. In order to simplify input, the cumulative loads 
(specified forces and displacements) at the end of each load increment 
are computed as the sum of the two load factors times their respective 
load vectors. Thus the two load vectors remain constant throughout the 
entire problem. For engine problems, the two vectors typically define 
the pressure distributions in the chamber and coolant tubes, respectively. 
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Although the present program does not allow the load distribution within 
a load vector to vary from one increment to the next, the two Itoad 
factors may of course vary independently of one another. 


C-4. Incremental Mechanical Load Data - These data are read by sub- 
routine READ4: 


NINCR 
NITER(I) 
PFACT(J,I) = 

CTIME(I) 


number of load increments to be run 

maximum number of iterations for increment I 

load factor to be applied to load reference vector 

J during Ith increment 

creep time increment for load increment I 


II-l. Thermal and Z-Di recti on Loads - These data are read by subroutine 
READ5: 

TIDENT = increment I.D. title 

T1(I) = cumulative (total) temperature at end of increment 

for element I 

ZSl(I) = cumulative Z load at end of increment for element I 

The increment title consists of any 80 characters used to describe the 
current load increment. Before the thermal and Z loads are read for an 
increment, T1 and ZSl are set equal to their values at the beginning of 
the increment (TO and ZSO, respectively). This provides default loads 
(unchanged from those of the previous increment) for those elements for 
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which the user does not specify a load value. The thermal and Z loads 
are then read according to data items I I- lb and c, in any desired element 
order. 

5.2 I/O AND STORAGE FILES 

The following are descriptions of files used in the BOP ACE program, 
listed by unit number. 

5 

Input card file. 

6 

Output printer file. 


UINI 

UINl is the user-defined unit number for type I input data. 

UIN2 

UIN2 is the user-defined unit number for type II input data. 

UOUT 


UOUT is the user-defined unit number for the major output data file. 


UINRS 

UINRS is the user-defined unit number for the input restart data file. 
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UOUTRS 


UOUTRS is the user-defined unit number for the output restart 
data file. 

UNITE! = 11 


File for storing the merged elastic stiffness matrix. 

UNITE2 = 12 

File for storing the decomposed elastic stiffness matrix. 

UNITP1 = 13 

File for storing the merged total Jacobian matrix. It is used only 
when the input variable SCODE is equal to 3, 4 or 5. 

UNITP2 = 14 


File for storing the decomposed total Jacobian matrix. It is used 
only when the input variable SCODE is equal to 3, 4 or 5. 

UNITS1 = 15 

Scratch file used for temporary storage by the Gauss wavefront merge 
and decomposition routines. 

UNIT52 = 16 

Same as UNITSl. 
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6.0 PROGRAM FLOW AND RESTART OPTIONS 

The major steps accomplished during a BOPACE run are shown in the program 
flow summary of Figure 6.0-1. 

Steps 3-7 accomplish the initialization of basic variables (program 
control constants, material data, mesh, and load vectors) and incremental 
variables, and the formation of merged and decomposed elastic stiffness 
matrices. These steps follow from reading of input data in the case where 
a new problem is being started, or from reading of the input restart tape 
in the case of a restart. If data is to be saved for a future restart, 
steps 9-10 write the basic variables and elastic merged and decomposed 
matrices onto the output restart tape. In step 11 the incremental 
mechanical load data, including load factors and creep time increment for 
each load increment to be run, are read. 

The remaining steps involve the incremental and iterative calculations. 
Updating of the Jacobian matrices in step 23 is performed only when 
convergence slows down, and in general this occurs only once per increment 
or once every several increments. 

In step 26 the computed incremental variables are written onto the output 
restart tape. If a future restart is provided for. This allows the user 
to request a restart from the end of any previously run load increment. 
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Figure 6.0-1: PROGRAM FLOW 
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Figure 6.0-1: PROGRAM FLOW (COMTINUFIi) 
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7.0 BOPACE ERROR MESSAGES 

BOPACE uses the FORTRAN STOP codes described in this section to indicate 
error conditions which may occur during execution of the program. Certain 
errors also generate a printed error message, in order to aid the user in 
locating the source of the error. 

Errors fall into two categories, those due to the problem definition or 
user input data, and those caused by a program or machine malfunction. 
Errors due to a machine malfunction rarely occur and in these cases a rerun 
of the problem will usually eliminate the error. If an error recurs and 
help is needed in correcting the problem, contact a BOPACE programmer for 
aid. Have available a listing of the input data, the printouts of the 
runs which failed, the input data deck, and a description of the problem. 

The following are explanations of the BOPACE error STOP codes, listed 
by subroutine in which they occur. 

READRS 

9999 Normal program exit caused by reading final blank card 
after all problems are run. 

101 Non-positive value input for a required file unit number. 
READl 

201 Plane stress-strain code input not equal to 0 or 1. 

202 Illegal value input for number of materials. 
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READTM 

301 Wrong material number input on material property card, 

302 Number of points input for material property curve 
exceeds maximum 

303 Non-positive value input for modulus of elasticity. 

304 Value input for Poisson's ratio is 0.49 or greater in 
plane- strain problem. 

305 No points input for a required material property curve. 
READTP 

401 Wrong material number input on plasticity type-code card. 

402 Illegal plasticity type input. 

403 Illegal kinematic hardening code input. 

404 Wrong material number input on plasticity temperature card. 

405 Plasticity temperatures not in ascending order. 

406 Number of temperatures for a material exceeds maximum. 

407 Hardening point defines non-positive yield-surface size, 
or negative kinematic value. 
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READTP 

408 Number of points input for a curve exceeds maximum. 

409 No points input for a required curve. 

410 First point input on a curve has non-zero abscissa. 

411 No curves input for a required hardening description 
of a material . 

READTC 

501 Wrong material number input on creep type card. 

502 Illegal creep type input. 

503 Number of points input for a creep reference curve exceeds 
maximum. 

504 No points input on the creep reference curve for a material. 

505 Wrong material number input on creep temperature card. 

506 Creep temperatures not in ascending order. 

507 Number of creep temperature factors for a material 
exceeds maximum. 

508 Number of creep stress factors at a temperature exceeds 
maximum. 
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READTC 

509 

510 
READC 

601 

READM 

701 

702 

703 

704 

705 

711 

712 

713 

714 

715 


No creep stress factors input at a temperature. 

No creep temperatures input for a material. 

I.D, of special coordinate system exceeds maxinfium. 

Mesh node I.D. exceeds maximum. 

I.D. of a node location coordinate system not equal 
to 0 or 1 . 

I.D. of a node displacement coordinate system exceeds 
maximum. 

Number of input nodes exceeds maximum. 

No nodes input. 

Element I.D. exceeds maximum. 

Illegal value input for element material number. 
Illegal node I.D. on an element. 

Number of input elements exceeds maximum. 

No elements input. 
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READ2 

801 Illegal node I.D. for a force- displacement-constraint 
specification. 

802 Illegal component number (not equal to 1 or 2) 

803 Illegal force-displacement-constraint code. 

804 Constraint specified with dependent (already constrained) 
DOF. 

805 Constraint specified between DOF on same element. 

READ3 

901 Number of input mechanical load reference curves not 
equal to 2. 

902 Illegal node I.D. on a mechanical load. 

903 Illegal component number (not equal to 1 or 2) on a 
mechanical load. 

904 Load input on dependent DOF constrained to DOF with 
specified displacement. 

READ4 

idol Number of load increments exceeds maximum per run. 
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READS 

1101 Illegal element I.D. on a thermal load, 

1102 Illegal element I.D. on a z-load. 

Linear Equation-Solver Routines 

5011 I/O Error during merge. Program or machine malfunc- 

tion during the merging of the elemental stiffness 
matrices to form the global stiffness matrix. 

5021 Bandwidth too large for decomposition save array. 

The bandwidth is too large for the amount of core 
storage allocated (see Equation 4.2-1). Corrective 
action: renumber the nodes to reduce the bandwidth. 

5023 No decomposition partitions available. Program or 
machine malfunction during decomposition. 

5024 Decomposition node not in active node array. Program 
or machine malfunction during decomposition. 

5031 Scratch array too small for solution work. Program 

or machine malfunction during forward and backward 
substi tution. 

5041 Illegal save tape I/O operation command. Program or 

machine malfunction in reading or writing the checkpoint 
tape. 
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7.0 . (Continued) 

Linear Equation-Sol ver Routines 

5042 Illegal matrix type. Program or machine malfunction 
in reading or writing the checkpoint, tape. 

5043 Illegal save tape defined for save operation. Program 
or machine malfunction in reading or writing the 
checkpoint tape. 

5051 Large decomposition not available. The bandwidth is too 

large to solve the problem using in-core decomposition. 
See Equation 4.2-2. Corrective action: Reduce the 

bandwidth by renumbering the nodes or reducing the 
problem size. 


7-7 



05 - 17266-2 


THIS PAGE LEFT BLANK INTENTIONALLY 


7-8 



PRECEDING PAGE BLANK NOT FILMED 


D5-17266-2 


8.0 REFERENCES 

1. R. J. Melosh and R. M. Bamford, "Efficient Solution of Load-Deflection 
Equations," Journal of the Structural Division, ASCE, April 1969. 

2. W. D. Whetstone, "Computer Analysis of Large Linear Frames," Journal 
of the Structural Division, ASCE, November 1969. 


8-1 



D5-1 7266-2 


THIS PAGE LEFT BLANK INTENTIONALLY 


8-2 



D5-17266-? 



9,0 EXAMPLE PROBLEMS 

Four example problems are provided here in order to demonstrate the 
BOPACE program capabilities and to aid the user in setting up his problems. 
The examples in Section 9.1 serve to demonstrate and check out various 
program options, including temperature-dependent elasticity, creep and 
prescribed normal loads. Section 9.2 illustrates a thermal ratchet 
behavior, which must be considered as a possible important effect during 
any combined mechanical loading and high- temperature thermal cycling. 
Section 9.3 illustrates the procedure for using actual cyclic test data to 
determine the isotropic and kinematic portions of hardening, and gives an 
indication of the test-analysis match which is possible with the BOPACE 
combined hardening approach. Section 9.4 gives an example of the required 
input data for a typical engine analysis, involving a relatively coarse 
finite-element mesh but with several different materials. 
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9,1 PROGRAM CHECKOUT PROBLEMS 

Plane Stress with Combined Hardening - The basic characteristics of 
BOPACE combined hardening are shown in Figure 9,1-1 for a uniaxial 
(special case of a plane-stress) problem. Figure 9.1-lA gives the 
assumed monoton ic stress-strain hardening curves. The size of the 
yield surface is defined by the isotropic stress* while the Bauschinger 
kinematic hardening is defined by the difference between the total 
stress and the isotropic stress. Thus the hardening is completely 
kinematic out to a strain value of 7.0 (plastic strain of 3,0), after 
which there are equal amounts of isotropic and kinematic hardening. 

A resulting cy cl ic stress-strain curve is given in Figure 9.1 -IB. The 
15 load increments were chosen so as to result in the exact a-e points 
given in the figure insert table. Note that the hardening parameters 
(k and K*^) in this example were based on effective plastic strain rather 
than on plastic work, because it makes the relationship between the 
monotonic and cyclic curves more readily apparent. 
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(B) STRESS-STRAIN PATH UNDER LOADING 


FIGURE 9.1-1: UNIAXIAL TEST PROBLEM FOR CYCLIC COMBINED HARDENING 
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9.1 (Continued) 

Plane-Strain with Additional Options - The plane-stress problem just 
described was altered to illustrate the use of several additional BOPACE 
options, including temperature-dependent elasticity, creep and prescribed 
normal loads. 

A BOPACE plane-strain checkout analysis was performed using the finite- 
element mesh and loading given in Figure 9.1-2. A listing of the input 
data and the printed output results are included at the end of this 
section. 


Y 



FIGURE 9.1-2:' PLANE-STRAIN CHECKOUT PROBLEM 
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9.1 (Continued) 

A summary of the problem is provided by Table 9.1-1. The 15 increments 
correspond to those of the previous plane-stress problem. The values 
of incremental plastic strain, stress, effective stress center, and yield- 
surface size given in columns 2-5 of Table 9.1-1 were kept the same as 
those of the plane-stress problem. The stress is equal to the product 
of the temperature- dependent elastic modulus (column 6) and the elastic 
strain (column 7). 

The creep strain listed in column 9 results from the material creep 
definition of Figure 9.1-3. There the reference creep curve for a 
strain-hardening material is assumed as shown in (A), while (B) defines 
the creep factor as a function of average, stress level during the 
increment. The creep strain may be determined using the average stress 
level (column 10 of Table 9.1-1), the creep factor (column 11) and the 
specified creep time increment (column 12). 

In addition, the Z-load strains given in column 13 and thermal strains 
in column 14 were imposed. In order to keep the results simple and 
exact (all numbers in Table 9.1-1 are given exactly), the Z-load and 
thermal-strain values were selected so as to give zero normal stress 
in each increment. For example, in increment 11 we have: 
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9.1 


(Continued) 
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-0.3 

- 1.0 

-0.5 

1.5 

-0.3 


Because the imposed Z-load strain also equals -0.3, a zero value results 
for the normal stress 0 ^ 2 * noted that this example can 

be used for either a plane-stress or a plane-strain checkout run. 


The prescribed displacements shown in column 15 were determined from the 
various components of the total strain. For example, again in increment 
11 : 
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FIGURE 9.1-3: CREEP DEFINITION FOR PLANE-STRAIN CHECKOUT PROBLEM 
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TABLE 9.1-1: RESULTS FOR PLANE STRAIN WITH Z-LOADS AND CREEP 


0 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

(S) 

( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

( 13 ) 

( 14 ) 

( 15 ) 

06 ) 

Incr . 

— 

°YY 

3 

f^YY 

|5-^i 


e 

"YY 

,P 

■'YY 

1 

"ave 



z - 1 oad 
strai r 

t 

c 

Oyy 

Temp . 

0 

- 

0 

0 

2.0 

1 .0 

0 

0 

0 


- 

- 

0 

0 

0 

1 .0 

1 

0 

2.0 

0 

2.0 

2.0 

1 .0 

‘0 

0 

1 .0 

1.0 

0 

0 

0.3 

1.3 

2.0 

2 

0.5 

2.5 

0.5 

2.0 

2.5 

1 .0 

0.5 

0.5 

2.25 

1 .0 

5.0 

0 

0.8 

2.8 

3,0 

5 

0.5 

3.0 

1.0 

2.0 

3.0 

1.0 

1 .0 

0.5 

2.75 

1.0 

0 

0 

1 .05 

3.55 

4 . 0 - 

4 

1.0 

3.5 

1.5 

2.0 

3.5 

1.0 

2.0 

1.0 

3.25 

1.25 

4.0 

0 

1.8 

5.8 

5.0 

5 

1 .0 

4.0 

2.0 

2.0 

4.0 

1 .0 

3.0 

1 .0 

3.75 

1.75 

0 

0 

2.3 

7.3 

6.0 

6 

0 

2.0 

2.0 

2.0 

2.0 

1.0 

3.0 

1.0 

3.0 

1.0 

0 

0 

2.3 

7.3 

7.0 

7 

0 

0 

2.0 

2.0 

1 .0 

0 

3.0 

1 .0 

1 .0 

1 .0 

0 

0 

2.0 

6,0 

5.4 

8 

- 0.5 

- 0.625 

1.5 

2.125 

1 .25 

- 0.5 

2.5 

1 .0 

0.3125 

1 .0 

0 

- 0.6 

1 .0 

4.0 

8.0 

9 

- 0.5 

- 1 .25 

1.0 

2.25 

1 .25 

-1 .0 

2.0 

0 

0.9375 

1.0 

10.0 

- 0.2 

0.5 

1 .5 

9.0 

10 

0 

1 .0 

1 .0 

2.25 

2.0 

0.5 

2.0 

0 

0.125 

1 .0 

0 

0 

1.15 

3.65 

10.0 

n 

0 

3.25 

1.0 

2.25 

3.25 

1 .0 

2.0 

1 .0 

2.125 

1 .0 

10.0 

- 0.3 . 

1 .5 

5.5 

4.5 

12 . 

0.5 

3.875 

1.5 

2.375 

3.875 

1 .0 

2.5 

T.O 

3,5625 

1 .5625 

0 

0 

2,05 

6.55 

12.0 

13 

0.5 

4.5 

2.0 

2.5 

4.5 

1.0 

3.0 

1.0 

4.1875 

2.1875 

0 

0.2 

2.5 

7.5 

13,0 

14 

1 .0 

5.25 ■ 

2.5 

2.75 

5.25 

1 .0 

4.0 

1,0 

4.875 

2.875 

0 

0 

2.8 

8.8 

14,0 

15 

1 .0 

6.0 

3.0 

3.0 

3.0 

2.0 

5,0 

4,125 

5.625 

3.625 

10.0 

0 

5.1625 

16.2875 

15.0 































16.297% 


10.0 


INCREMENT 

i 




10 

2.0 

, ?0 2.0 



increment 

^ 




10 

0 

20 .3,0 



INCREMENT . 

.. . _3 




1C 

4,0 

20 4.0 



INCREMENT _ 

4 




10 

5.0 

20 5,0 



INCREMENT 

5 




1C 

6. C 

20 6.0 



INCRfl4£ftlI__ 





10 

T.C 

20 7,0 




7 




if 10 

> 

5.4 

20 5.4 


•K4 

ro 


p 



cn 

t 

fO 

1 c 

R. 0 

20 R.O 



10 

-0.6 

20 -0,6 



increment 
L£L 



20 ^.0 



10 

-^L.^ J 

_20_- Q. 2 



XKCRfMBlT . 

.. .. 1.0. 




10 

10,0 

20 IC.O 



10 

0. 0 

2C 0.0 



INCREMENT ■ 

IQ_ 

11 

-70 4. .5 ■ .... 



iSL 

-0.-' 

20 -0..3. 



lACREMENT., 





10 

1 2, 0 

20 12.0 



10 

0.0 

20 0 .0 



INCREMENT 

LO. 

13 

13,0 

20 13.0 




Q_ U tO. 



10 0.2 


0.2 


I.MC«EMENT 

10 ] i.O 


1 

:^o 14.0 


ic: 0. c 


0 . 0 


INCPFKfrMT 

10 . 15,^1 


I 

t 



I 



?C 1^.0 



U-1’6 


SXaPTlNC PPCBLEK 


.0 C P.AC E _G ej^.P 4 L C HECK {Pl4^E STRAIN »ltTH Z~LCACS« EL4ST IC-PL AST I C-CR f EP ) 06/08/73 

SClUTIC^ f'FTHCC CORE =_ 5_, 

^AXIMUP KC. STIFFNESS UP’CftTES PER INCREMENT » ' 3 

MAXIMUM TDTAL ITFPATJ.CNS PEP INCRE MENT = 10 _ 

MAXIMUM ELASTIC ITEftATICNS PEP INCREMENT = 2 

MAXIMUM MAGMTUCE FCP ELASTIC-PLASTIC SUM CODE « 2 

MAXIMUM RECCCTICNS » ' 1 ^ ^ 

CCNVERGENCE RECUCTICN FACTCP « ' O-SCOCOF 00 
maximum' SPECIFIEO EPBCP'NGPM s O.IOOOOf-OA 

fPACJK.N FRCM end of INCREMENT TO EVALUATE SLCPE « O.IOOCOE 00 

PLANE-STRAIN PROBLEM 

NO, OF materials'"^ 1 “ 

DEFAULT TMC_K^S => C.IOOCOE 01 

FABPICAtiCN TEMFERATUPE s C.IOOOOE 01 


MAT£FJ A.k_NC_. _ 

1 TEMFERATURF CEPENCENT PROPERTIES 





O.IOCCE 01 
e.20flCF 01 

0.0 

0-3CC0E 00 


C.3000E Cl 
CiAOOOC U 

o.accoE 00 
0.1050E Cl 


C.A500E Cl 
«> C.500CE Cl 

0.1500E 01 

o_.Lec_oj_oi 

o 

trr 

C.60CCE 01 
h. C.700CE Cl 

C.23C0E 01 

•VI 

J\> 

2? 

C.8Q0QE 01 
C.9000E 01 

O.IOCOE 01 
0.50C0E OC 

» 

ro 

C.lbOQE C2 
C.12CCE C2 

0.U50F 01 
0.2C50E 01 


0.130QE Q2 
C.IACOF C2 

0.2500E 01 
0.28CCE 01 


C.ISCCE 02 

0.S162E 01 


temperature 
C.lflCOE Ql 

ELASTIC MDC. 
O-IOCCF 01 


C,2000F 01 
C,3000E Cl 

0.20CCF 01 
O.Z5COE Cl 


C.AflCCE Cl 
C-5000E 01 

0.3CC0E 01 
0.35C0E 01 


C.5A0CE 01 
L.6_0CCE 0.1 . 

O.ICCCE 01 
O.ACCOE 01 


C,7QCCE Cl 
C,.800CE .01 

C.2CCCE 01 
0.1250E Cl 


C*9Q0CF Cl 
.O.IOOOE C2 

0.I250E 01 
0.20COE 01 


C.12CCF 02 
C.1300F 02 

0.3675E Cl 
0.A5OOE 01 


C.1400E 02 
. C,J50OE. 0,2 .. 

0.5250E Cl 
0.30CCE 01 


temperature 
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b.abooF 00 
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0.350CCP 01 


PAPAMETEP ■ BIKEBATIC SHAPE 

O.C 0.0 _ 
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0.0 o.o' 

_C_.l_QOCpE Cl C.IOOCCE 01 
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hC. Cf, LC^C INC_BEHENTS_ = _ 15 

INCREMFKT MAX, fTERATICNS ^ECHAMCAL CURVE FACtCRS CREEP TIRE 
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INCPEMEM 




CUMULATIVE thermal LCACS 

AND Z-LCACS FCR LCAC INCREMENT 1 

INCFEMEM 


1 



ELEWEKT 1. 
TEMP. 

C. 20 

C.20000E 

10 

Cl C.20QCOE 01 . 



Z-LCAO 

0,0 

c.o 



ITERATIVE 
ITERATIVE 

ERRCP » 
ERRCP s 

0.7A957E 00 
O.S1221E 00 



ITERATIVE 

ITERATIVE 

ERRCP = 
ERROR = 

0.7A957E 00 
0.2C911E-Qfc 



I tepative 

EPRCR = 

0.18C29E-06 















END 

C F L C 

AC INCREMENT 1 

INCPEMEKT 


1 



MECFANICAL 
C(^P TIME 

LOAD CURVE FACTORS * C, 

INCREMENT » 0-0 

13C0E Clt 

0.0 

NC. ELASTIC ELEMENTS * It NO, PLASTIC ELEMENTS » 1 

1 ELEMENTS HAVE CHANGED ELASTIC TO PLASTIC, 0 ELEMENTS PLASTIC TC ELASTIC DUPING THIS INCREMENT 

SPECIFIEC 
SPECIFICC . 

MAX, NO, 
MAX. NO- 

STIFFNESS UPDATES = 3, 

ITERATIONS PER UPDATE = 

NC, UPDATES PERFGPMEC * 0 

1C. NC. ITERATIONS PERFGPMEC SINCE LAST UPDATE * 5 


specified wax, UNeALA^CEC•-FCPC€ ERRCR = O.lOCOE-04, ACTUAL ERROR * 0.iB03E-06 

L__. *»»** »» C U PUIATIVE INTERNAL FORCES AND C I SPL ACEWENTS ******** 

NODE *• ♦*•**♦*♦* forces **♦**•*•♦ ***** CISPLACEHENTS ****** 

NO. _i , c , u y _y y 

... _ l 30 0.23a4lE6 E-C6 0.99S9S9CE 00 0. 172L91 LE«Q6 C,12999^q E 01 

2 4fl -0.I976'633F-0« C.IOCOCCCE Cl 0. 38«a754E-C6 0.1299999E 01 ' 

3 1C -0,IU62A^-Q_6 -0,999999CF 00 ^0 OiO 

A 20 6,73869C5E-C7 -C.IOCOCCCE 01 ' 0. 5a8856CF-06 O.C 
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9.1-16 


»»*»* »»» ♦»i>» **» »**»»»»** * *»**»»*»» ** _6LA^TIC STRAINS *♦»*»«♦»*»♦ *»* *♦♦**♦*♦»**»**< **» ***» *»* 

ELfMEKT ♦** THFPt^Al STRAINS ♦**********i** INCREMENTAL ♦*♦♦*****■**•***•* ************* CUMULAfiVE ♦*♦•*****•*#**♦•** 


NC. 

I .C. 

INCREMENTAL 

CUMULATIVE 

XX 


YY 

ll 


XY 

XX 


YY 

ll 


XV 

1 

20 

0.300CE CC 

0. 3000F 00 . 

C.IOOOE 

Cl 

-C.3000E 

CO -0.3000E 

00 

0.1020E-06 

C.IOOOE 

01 

-0.3000E 

00 -0.3000E 

00 

0.1020E-06 

2 

10 

0.3000E CC 

C.30C0E OC 

C.3500E 

CO 

C.3500F 

CO -0.3000E 

00 

0^650CE 00 

0.3500E 

CO 

0.3900E 

00 -C.30CCE 

00 

0.65C0E 00 


^ »*♦» »»♦»♦»»»*»♦♦»»»<■»»»»»»»»♦ plastic strains ******************************** ♦»»»»»«> 

ELEMENT ’ ***'**"PLASTIC WCPK **** »*»•'****♦••** 'INCREMENTAL ♦***♦♦•***♦*♦**♦♦ CUMULATIVE ♦♦**••♦♦*•**♦♦♦♦♦* 


KC. 
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CUMULATIVE 

XX 

YY 

ll 

XY 

XX 

YY 

ll 

XY 

1 

2C O.C 

0.0 

C.O 

C.C 

O.Q 

0.0 

C.O 

0.0 

C.O 

0.0 

2 

ic 0.0 

c.c 

c.c 

C.O 

0.0 

0.0 

C.O 

0.0 

O.C 

~o.b 




1 20 Q.G__ 0. 2000E □ I 0*0 C.O 0.0 0,0 C«2000E 01 -0»9 170 E-Ca ~0»^ 677E~C6 0«1S 6^E ~06 

2 IC O.O” 'C.2000E Cl C.O 0«0 0,0 0.0 C.IOOOE 01 C.IOOOE 01 -C.68‘77E-C6 C.IOOOE 01 


ELEMENT E-P SUM I NCREMENT AL JXTAL SUREACE *»*» EFFECTIVE PLASTIC STRAINS **» ***** EFFECTIVE CREEP STRAI NS ♦♦»* 

NO. I. D. CODE CODE TEMPERATURE t‘emPE“rATUPF ' Y 1 ELC SUE INCREMENTAL SUM INCR. CUMULATIVE INCREMENTAL SUM INCB. CUMULATIVE 

1 20 1 1 C*10C0E 01 0.2CC0E Cl C.2CCCE 01 O.Q 0^0 0^0 0_t0 0^0 O^C 

2 ~ ic 0 -2 O.IOCCE 01 C.2CC0E 01 0.200CE 01 0,0 C.O 0.0 0.0 C.O O.C | 
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CLMLLATIVE THERMAL 

LCACS ANC Z-LOACS 

FOR LCAC 

INCREMENT 2 


INCREMENT 


2 





ELEMENT !, 
TEMP. 

.C. 20 

C.30000E 

ID 

Cl C.30000E 01 





2-LCAD 

c.o 

C.O 





ITERATIVE 

ITERATIVE 

epr'cr » 

EPfiCR « 

0.S2167E 00 
o.a^eiAE 00 





ITERATIVE 

ITERATIVE 

EPRCR = 
ERRCP = 

O.esSlAE 00 
O.Se787E 00 





ITERATIVE 

ITERATIVE 

EPRCR = 
EPPCR * 

0.EU47E CG 
0.E8986E 00 





ITERATIVE 

ITERATIVE 

ERRCR = 
EPRCR = 

O.S8201E 00 
0.61647E CO 





ITERATIVE 

ITERATIVE 

EPRCR = 
ERROR » 

0.e8986E 00 
0.9E201E 00 





ITERATIVE 

ITERATIVE 

EPRCR = 
ERROR = 

0.E1647E 00 
0.E5C69E 00 





ITERATIVE 

ITERATIVE 

ERROR = 

error = 

0.76A68E OC 
O.AA035E 00 





ITERATIVE 

ITERATIVE 

ERROR » 
ERROR « 

0.3e2lAE-01 

0.35867E-02 





ITERATIVE 

ITERATIVE 

ERRCR » 
ERRCR a 

0.3A599E-03 

O.33219E-0A 





I TERATIVE 
I TERATIVE 

EPRCR = 
ERROR = 

0.37867E-05 

0.92718E-06 











o' 

yi 







rv5 



END C 

f L C A C INC 

P E M E N 

T 2 

Ol 

t 

ro 


INCPEH^Nt ""2 

MECt-ANKiL LOAD CURVE PACTORS = O.Z8COE Oli 0.0 
CREEP TIPE increment = 0.5000E 01 

NC. ELASTIC ELEMENTS = 0, NC. PLASTIC ELEMENTS » 2 

1 ELEMENTS HAVE CHANGED ELASTIC TC PLASTIC. 0 ELEMENTS PLASTIC TC ELASTIC DUPING THIS INCREMENT 
SPECIFIED MAX. NO. STIFFNESS UPDATES = 3. NC. UPDATES PERFCPMEC » 1 

SPECIFIEC MAX, NO. itERATICNS PER UPCATE = 10. NC. ITEPATICNS PEPFGRHEC SINCE LAST UPOME * 10 

SPECJFIEC MAX. UNBALANCEC‘ F0PC E ERRTR = C.IOC OE -OA, A CTU AL ER ROR = 0.92 72 E‘-0 6 


■>*»***♦ CUMULATIVE I NTERNAL FO RCES ANC C I S FL ACF M ENT S *»■»»*** * 

•• NODE ** •*••*•**♦ FCPCES ♦*•*♦*»** »*♦** DISPLACEMENTS *♦»*** 

- U.Ct u y y 


1 30 “0.A172325E-C6 0. 12AqS9SE 01 C . 206 l fl2 3E-C6 C.2799SS9E Cl 

2 AO 0.72i9A3aE-07 0, 12<»9999F 01 0. 7573697F-C6 0.2799999F 01 

3 1C 0.S1122A1E-06 -Q. 12A9S99F Cl 0^0 C.C __ _ 

A 2C -O.U6lfl60E-06 -0.1250CCCE 01 C.89A5562E-06 C.C 



9.1-18 


***♦*»<< *»»»»»» »»»»*♦»»» ELAST IC STRAINS ♦**♦»»»**»♦♦*♦ **»»********»******i»**«»i> 

ELEMENT *** THEPM4l STPAINS **♦ ****j»»******* INCPEMENTai •♦•••♦*•♦*♦****•* "'«»*i*****»4<**4< CUMULATIVE’ 

M . I .C .„I NCPEMENTA L CUMULATIVE XX VV 11 KV )0< 11 XV 

1 2a O.SCOaE OO C.eOOCE OC • C >77A9E*06 -C. B345E-C6 Q»41?2E-06 Q .a8 9»E-0B O.IOOOE Cl -0.30C0 E 00 -C.3 0CCE CO 0-1109E-06 

2 10 0.50COE CC C.eOCOE 00 -C,1780E“C6 C.4768f-C6 O.C Ci23aAE-06 d^35QQP CO” 0.35C0E 00 -C.iOCOE OO’ 0.6566 e'60 


. »»»»»»»!»» ♦♦» »»»»»♦»»»»♦<»#»» »»»»»»»» PLASTIC STRAINS •♦****••*♦**♦*♦♦ t^***!!!***************** 

element PLASTIC WORK ♦******♦♦♦*♦* INCREMENTAL *•#*•**•**•♦**<*** CUMULATIVE ** * aT'^ ^ 

iLLt. _-l -J* -LN CREHF NTAI CUMULATIVE XX YY ll XY XX H U *1. 


1 20 0.U2 5E Cl_ 0«.1125E 01 C.SOOOE CO -C .2500E CO >C.250CE 00 -0. 256S E-07 C «500Q€ CO -0.2 5C0E CO -C.2500E 00 -O.2 505E-O7 

2 10 0.U25E 01 0,11256 01 0.12506 CO C.1250E CO -0.2500E 00 0.3750E CO 0.1250E 00 0.1250E 00 -0.2500E 00 0,37501" 


****» **» »»♦♦»»»*♦*»»»»**♦♦♦»*»*»*»♦» C REEP STR AINS 

element •*♦*#* CREEP WORK *•»*• ************* INCREMENTAL *•*****♦***♦»•*♦* «********«**4r CUMULA'tIVE 

NO. I .C. I NC P EM ENTAL CUMULATIVE XX YY ll XY XX YY Li XY 

1 20 0.1 1256 01. 0.1125F 01 O .SCOCE CO -C,2500E CO -0.25QC E 00 -0 . 164AE -07 C. SOOOE 00 -0.2500E OO -0.25006 00 -0,16««E-07 

2 10 0.11256 Cl C,1125E 01 C, 12506 00 G. 12506 00 -0.25006 00 0.3750E 00 0.1250E 00 0.1250E 00 -C.2500E OO O,3750E 00 


CUMULAT IVE STRESS QU ANTITIES ♦♦**•♦*****!***•**•#*♦♦**♦♦♦♦♦**♦* 

ELEMENT EFFECTIVE EFFECTIVE ****♦**«►*'♦** STRESS CENTER **************** *************** STRESS ♦•♦•••********^*#** 

NO. LiIj^ C£MfP STRESS XX YY Zl XY XX YY ll XY 

1 2Q 0.5Q^E CC_ .Q,250CE__0.1 0.3333E CO -C.1667E QO -C.1 667E O C -0.1723E-Q7 C.2500E 01 -0.5456E-06 0.132CE-C5 C. 21326-06 

2 10 0.50006 CC 0.2500E 01 C. 03336-01 C.8333E-01 -0.16676 00 0.2500F 00 C. 12506 01 0.1Z50E 01 -0.12846-05 O'iTISdT'Cl 


ELEMENT E-P. .SUM INCR EMEN TAL lEtAL SURFACE »»»♦ EFFEC T IVE PLASTIC STR AINS *** * **** EFFECTIVE CREEP STRAINS 

NC. I.C. CQCE CGDE TEMPERATURE TEMPERATURE YIELD SIZE INCREMENTAL SUM INCR. CUMULATIVE INCREMENTAL SUM INCP. CUMULATIVE"' 

I — ,.|C 0 2 _C.10C0F 01 Q.30QOE 01 Q.200QE 01 Q.5000F CO 0.50QC6 00 0. SCOPE 00 0. 5QOOE 00 C.50C06 CO C.5QC06 QO 

2 10 I 2 O-IOCOF 01 C,3CC06 01 0.20C0E 01 O.SOCOE CO 0,50006 00 0.5000E 00 0.500QE 00 C.50CC6 00 0.50006 CO 


D5-1 7266-2 





6l-L'6 


3 


IKCPEMi’Kf' 3 

lo To 

_ .. TEMP. C.4.Q0.CLe.E Cl C.40000E 01 

Z-tCAD C.O C.O 


ITEPATIVE EPBCP = 0.eil6CE 00 


ITERATIVE ERPCP 


ITERATIVE EPRCP = 


ITERATIVE 6PPCP 


s 

0.SZ446E CO 
C.34C5RE 00 




1 

1 

* 

0,lfc222E CO 
C.46C23E-01 




- I 

s 

0. 16448 E-Ol 
0.52C84E-C2 




■ 


0.17755E-C2 

0.56S27E-03 






0, 28329 E-C4 
0.2C166E-05 





= 

0.U412E-05 






E N C C F LCAC INCREMENT 3 


JiiCREMEM_ _ 3 _T 

MECMANICAI LOAD CURVE FACTORS = 0.3550E 01, 0.0 ^ 

CREEP T IME INCREMENT » 0.0 ; m 

NO. ELASTIC' elements" « 0, NO, PLASTIC ELEMENTS -2 T 

0 ELEMENTS HAVE CHANGED ELASTIC TO PLASTIC, 0 ELEMENTS PLASTIC TC ELASTIC OURING THIS INCREMENT 

SPECIFIEC MAX. NO. STIFFNESS UPDATES = 3, NC. UPDATES PERFCRMEC = 1 

SPECIFIED MM. NOi ITERATICNS PER UPDATE = 10, NO, ITERATIONS PEPFQPME C S INCE LAST UPDA TE = 3 

SPECiPlED max'. UNBALANCE C-FOPCE ERRCR = C.lCCOE-04, ACTUAL ERROR « 0,174iE-05 


• ** CUMULATIVE INTERNAL FORCES AND C 1 SPLACEM ENT S **»»»»»» 

•• NODE ** **•**♦*•» FORCES *4**A***A ***•• DISPLACEMENTS ****** 


NO, 

I.D. 

u 

V 


u 

V 


1 , 

2 

30 

0.596C470E-C7 

0-1499999E 

01 

0. 1926832E-06 

C. 3549999E 

Cl 

40 

-0.4837949E-06 

0.1499999E 

01 

-0.92045I9E-C6 

0.3S49999E 

01 

3 

1C 

-0.263635EC-C6_ 

-0.1499999F 

01 

C.C 

C.C 


4 

20 

C.6878261E-06 

-O.I5COCOCF 

01 

-0.6157362E-C6 

C.C 



I 






02-re 


ELASTIC S TRAINS ♦ ♦♦♦«»» »*** **» *»*»»»»»***»*i» »**»»♦♦»»*# 

ELEMFM •** THEPM4L STPAIN$~*»* ************* IKCPEf'ENTfll *♦•*♦******»*•**♦ ' ******** CUMULATIVE ♦♦*•***•*♦••♦♦**♦* 


NC. 

I .C. 

incremental 

CUMULATIVE 

XX 

YY 

II 

XV 

XX 

YY 

U 


XY 

1 

2C 

0.25CCE CC 

0.1050E 01 

C.5960E-C6 

-0.2384E-C6 

0.2980E-C6 

0.5012E-07 

C.IOOOE 

01 -0.30CCE 

CO -C.300C6 

00 

0. 16106-06 

2 

1C 

0,2500r CC 

0.1050F 01 

-G,2384E-C6 

-0, 178BE-C6 

0,83456-06“* 

0*7I53E-06 

0.3500E 

00 0.3500E 

OO -0.3000B 

00 

0.65006 00 



ELEMEKT PLASTIC ViOPK *♦** 

>0. l.C. INCPEMgNTA L CUMULATIVE 


* *********************************** plastic s trains ♦o*******i*************i***********%**«** 
****♦**»•*♦•* INCPEMENTAL ♦***♦•*♦*•**♦♦**♦ 4********«*** CUMULATTvF •*****’iM*ir***' 

YV IZ XY )LX YY 1 1 XY 


2C 0^13T5E Cl 0.2500E 01 

10 0.1375F Ci C.2500E 01 


O.SQOOE CO - C«2 5 00E CO -0«2500E 00 0.3 365 Ej^OT. 
C.1250E CO 0.125CE 00 -Q.2500E 00 0.3750C 00 


C«1000E 01 -0.S0 CQ6 00 -0.5000F 00 0,77?9|r08 
C.2500E CO 0.2500E OO -0-5000E 00 0.75006 00 


******* ** ******************** CUMULATIVE STRES S_CU^N1XT lES 


ELEMENT 

EFFECTIVE 

EFFECTIVE 

STRESS CENTER **♦***•♦•*•***•* 

*4) *««***« STRESS 

******************** 

NO. 

I.C. 

CENTER 

STRESS 


XX 

YY Zl 


XY 

XX 


YY 

Zl 

XY 

1 

zc 

O.IOOOE OL 

C.3000E 

01 

0.6667E 

CO -0.3333E CC -0.3333E 

OC 

C.5199E-08 

C.3000E 

01 

0.1004E-05 

0.4420E-05 

1 

0,37156-06 n 

2 

10 

O.IOOOE 01 

C.3000E 

01 

0.1667E 

00 0.1667F 00 -0,33336 

OC 

0.5000E 00 

C. 1500E 

01 

0.15C0E 01 

0.6236E-06 

0.15006 Cl ‘f 1 


ELEMBAT L“P _ SUM INCPEMFNT^ _ TCTAL SURFACE *♦** EFFECTIVE PLASTIC STRAINS *♦* ***** EFFECTIVE CREEP STRAINS 

NO. I.C, COCF CODE TEMPERATURE TEMPERATURE YIELD SIZE INCREMENTAL SUM INCP. CUMULATIVE INCREMENTAL SUM llvCR^ CUMULATIVE 

t 20 C 2 C.IQOOF 01 Q.AOQQE 01 0.2Q0QE 01 Q.5Q0QE CO 0. I0Q06 01 O.IOOOE Cl 0.0 C.50CQE 0 0 0.5000E 00 

10 0 2 C.IOOOE 01 O.AOOOE Cl O.ZOCOE 01 0.50CCE CO O.IOOOE 01 O.IOOOE 01 


D5-17266-2 





C U f<u L A T I VE_ iHf? Mi LOACS ANC Z-LCACS FC«» LCAC IKCPEyFAT f 

INCPEKEM A 


EtEWEAiT 

TEMP. 

Z-ICAO 


I.C. 20 10 

C.5000CE 01 C. 5CQC0E 01 

0,0 c.o 


ITEPATIVE eppcp 
ITERATIVE EPPr.P 
ITERATIVE ERRCR 
t TEPATI VE ERRCR 
ITERATIVE EPRCP 
ITERATIVE EPR£.R 
ITERATIVE ERRCR 
ITERATIVE ERRCR 
ITERATIVE ERRCR 
IIEMTIVt ERRCR 
ITERATIVE ERRCR 
ITE RATIVE EPR^P 
ITERATIVE ERRCR 


« O.SCCtlE 00 
° Q,^C27CE 00 
= 0,61C26£ 00 

* C,eCi.8C_E 00. 

* O.lklCE 00 


, O.EESqT E OQ 
* 0.3C772E CO 

= 0.£C7<t7E’Cl 

= O.HS'SOE-Cl 
=_ 0,27Ee6E>02 
» 0.«i28SE-C4 

Jf 0. E5623E-05 
= '0.tC714E-05 


J 


ENC CF LCAC INCREMENT A 


. -IMIRE.HE>T _ A 

MECFANICAL LCAC CURVE FACTCRS = O.SeOOE Oii 0.0 

^ CREEP TIME INCREMENT = C.4000E 01 

-• KC. ELASTIC ELEMENTS * 0, NO. PLASTIC ELEMENTS * 2 

0 ELEMENTS HAVE CHANGED ELASTIC TO PLASTIC. 0 ELEMENTS PLASTIC TC ELASTIC DURING THIS INCREMENT 

SPECIFIEC MAX. NO, STIFFNESS UPDATES * 3, NC. UPDATES PERFCPMEC * I 

tPfCJPIEO. MAX. NO, ITEPATICNS PEP UPDATE = I0i NO. ITERATIONS PEPFOPMEC SINCE LAST UPDATE ■= 3 

SPECIFiEC MAX. UNPlALANCED-FCRCE ERPCR = 0. lOCOE-04, ACTUAL ERROR = 0.6071E-05 


• * • * * ♦. * . .,C UMULATIVE INTERNAL FCRCES AND DISPLACEMENTS ******** 
♦ • NODE ** ********* FCPCES *♦**♦♦*** ***** OISPLAC'EMENTS ****** 

.-..NiOi_ I.C. D y y y 

1 3C C.1490116E-05. 0.174^96E_01 -0. 5040 76 6E-06 0.5 799SS9E Cl. 

'-0.246i763E-d5 0.i7’49997E 01 - 0 . 1485 203E-O5 0.57999S9E O'l 

C.63Q4343E~C6 -C.I 749995F OJ £iO OjO 

0.34121346-06 -0.1749999F 01 -0 . 79572 1 5E-C6 ( 


20 




\ 


fLEKEM 
AC.* 1 . ti 


♦** THEPMiL strains *** 
incremental cumulative 


t t*rn***** ********** *****»***'9******* ELASTIC ST RAIN S »»♦ *♦»*»*♦♦♦♦*♦♦♦»*»*» »»*»»*»»»♦»»♦*»♦♦ 
INCPSRENTAL *****************' ' ************* CUMULATIVE *•♦******♦♦****♦•• 

XX YY ll . XY XX YY U XY 


1 20 0.75C0ECC O.IBOOE 01. 0.2384E- C6 ~0«6557F-C6 0.2384E-06 0.2066E -06 C. lOOOE 01 -C. 3 0CC E C C ~ C.3000E 00 0. 3676E "C6 

2 10 0.75C0F CC 0.1800E 01 -0«29a0E-C6 -C,2SflOE-C6 -0.3576f~0fc 0.3576E-06 C.3506 e CO 0.3500E 00 -0.3000E 00 0.6500E 00 


I 



******************************* ***** plastic strains •**♦***♦*♦**♦♦♦***♦*»**•♦•••*•**•♦**••• 

ELEMEM PLASTIC MCRK •*** ************* INCREMENTAL ***************** ' ************* CUMULATIVE *•#♦♦*♦•*#•»♦••••* 

NO. J_. C. I NC R EMENTAL CUMULATIVE XX YY ll XV XX y/ ll XY 


i 20 Q.3250 E Cl 0.57SQE 01 C. LOOOE Cl -0.500QE CO ~0.5000E 00 -0.1982E-07 C.20C0E 01 -O.IOOQ E 01 >C.10QCE 1202E-^ 

? 10 0.3250E Cl 0.5750F 01 C.2500E 00 0.2500E 00 -0.5000E 00 0.7500? 00 C.5000? 00 0.50COE 00 -C.lbbOE 01 a.lSOOE 01 



_ _ CREEP strains »»»«»»»»»»»«♦♦♦♦♦»«♦<♦♦♦»♦♦*»♦»»»»♦♦»♦<» 

ELEMENT •*****~C^EEP WORK ***** ************* INCREMENTAL ♦**•*•♦**♦♦*♦•♦** ♦♦♦*******♦♦* CUMULATIVE ♦*•*•********•*♦•* 



***************************** CUMULAT IVE S T RES 5 QUA NT I TIES ♦»♦»*♦*»♦» »»♦»♦»♦♦ » ♦♦♦♦*«* »** ♦»*« 


ELEMENT 

EFFECTIVE 

EFFECTIVE 

1A*»»**»*»**« stress 

center **************** 

**»*«*****•**** STRESS »»*«***MM*i*****4«««« 

NO. l.C. 

CENTER 

STRESS 


XX 

YY 

ll 


XY 

XX 


YY 12 

XY 

1 20 

0.1500E 01 

0.3500E 

Q1 

G. lOOOE 

Cl -C.SOCCE 

CO -0.5COOE 

00 

-0.14CeE-08 

C.3S00E 

01 

-0.3072E-C6 O.A062E-05 

0.9896E-06 

2 10 

C.1500F 01 

0.3500F 

01 

0.2500E 

00 0.2500E 

00 -0.5000E 

00 

0./6C0E 00 

G.1750E 

01 

0. 1750E 01 -C, 108CE-C5 

0.1750E 01 



F - p ___sy M I NCREMENTAL T C_TA L SURF ACE ^*_£ Ff_E CU Y3_ ^ A SJ_I C_ S IF A J_N S_* ** ♦***♦ EFFECTIVE CREEP STPA IN S **♦■«■ 

CODE CCDE temperature TEMPERATURE YIELD SIZE INCREMENTAL SUM INCP- CUMULATIVE INCREMENTAL SUM INCR. CUMULATIVE 

0 2 G. lOOOE 01 0.50COE 01 0.2000? 01 .0. lOOOE 01 0.2000E 01 0-2000E 01 C.SOCCE 00 C.IOCOE Cl C.IOCOE 01 









9.1-23 




CUKLtA TtVE THEPf^AL LOaCS />NC Z-LCaCS fC P LCA C ^^ CPE>'E^T 5 

5 


eleVekt i.c. 20 10 

TEMP. C.60000E Cl C.6000Q E 0 1 . 

l-LCAO C.O C.O 


ITERATIVE EFPCR = 0.65A54E 00 
ITEPATIVE EPRCP » 0.tlS6'%E 00 
ITFPATIVP FPRCB = 0.S3655F CO 
1 TEBATIVE EPPCB.i... 0.632 8aE 00 
iterative EPPPP » 0,2910«E 00 


..ITEPATIVE. EPRCR a 0 .382S9E~Q1 
ITERATIVE EPPCjR * d,6C3S6E-02 

ITERATI V E EPRCR « O.S2193E-03 

ITERATIVE FPRCP = O.UlOOE-03 
.ITERATIVE E RROR ^ 0.2I530E‘04 

fTEBATIVE EPRCP = 0,457i3E“06 

. I..llPA.Tlyf ERROR « 0.34204F-05 


ENC CF LCAC 1^CREM E_^LJL 5 


1NCPEMEA7 5 ' " 

...PEC.MN I CAL LOAC CU RVE EACTCP S - 0.73CQE 0 1 , 0.0 

CREEP TIME IMCPEMENT = 0.0 

_.lyie.... J LAST K... ELEMENTS =._ C. N O. P L ASTIC ELEMENTS = 2 ^ 

0 ELEMENTS HAVE fHANGEC ELASTIC TC PLASTIC, 0 ELEMENTS PLASTIC TC ELASTIC DURING THIS INCREMENT 

SPEClFlfC MAX. NO. STIFFNESS UPDATES = 3, NC. UPDATES PERFORMEC = I 

SPECIFlEf MAX. NO. ITERATIONS PEP UPDATE = 10, NO. ITERATIONS PEPFCPMEC SINCE LAST UPDATE » 2 

SPECIFIEC MAX. UNBALANCE C-FCRCE EPRCR = D.IOCOE-OA. ACTUAL E RROR = 0.3A2QE-05 


**** *** CUMULATIVE INTERNAL FCRCES AND C 1 5PL AC E MENTS **»■»*»*» 

NODE •* ********* FCRCES n******** ***** DISPLACEMENTS ****** 

NO. I.C. y y y y 

1 30 -0.4172325E-06 0.19SSS95E Cl ~0. 17823C9E-05 0.7299 9 <)9E 01 

2 40 -0.e503912E-C6 C.1999996E 01 -0. 196t750F- C5 C.7299999E c'l 

3 1C 0j_24ia51OE-C5 -0. 1 9999 94E 01 0^0 0^0 

4 20 -0.il‘5094tE-C5 -0,l99999eE 01 -C. 123E942E-C5 O.C 



9.7-24 i 


*■$*******************************■* ** El AST 1C STROKS »♦* ♦♦»»** » »**♦*♦ <»» »»♦♦» »»*»**»»^****<i» 

FLEMEM *** THERMAL STRAINS •** ♦**♦♦*#*****« INCREMENTAL ♦***♦**♦♦***♦****’ ' ■♦♦♦**■*••*'♦*** CUMULATIVE **»**»'***** **#*•*♦ ' 


-AC* L &£.±... 

INCREMENTAL 

CUMULAT I VE 

XX 

YY 

AL _ 

XY 

XX 


YY 

U 

XY 

1 20 

0.5000F CC 

ai.23Jl£._0L;_ 

0.17B8E-G6 

“0.9S37E-D6 

a.35?eEr0fc 

C. 1464E-06 

O.IOOOE 

01 

'0. 3 CCOE 

CO -C.30CCE CO 

0.5140E-06 

2 1C 

C.5000F OC 

0.2300F 01 

C.8345E-C6 

0.7153E-C6 -O.0345E-C6 

0,59606-07 

C.3500E 

cb 

0.35CCE 

CO -C.30CCE CO' 

C.' 6 5 COE do 


♦♦»»»»»»»♦»♦»♦♦»♦»♦»»♦»♦»♦»♦**»«»»*♦ PLASTIC STRAINS _ ♦* * »»» ? **♦**♦ *»» » ** »»*« ^ 

ELEMENT **♦♦♦ Pt"ASTIC WORK #*#♦*•♦♦****♦ INCREMENTAL •**»*♦**♦♦*** CUMULATIVE ♦♦♦'*♦•*'* ***♦'*•* ” 

NC. uc. incr emental cumulative %% >v u XV XX yy Zl XY 

1 20 0,375<3E CL OiSSOOE C,10 00E 01 -O.SOOOE 00 -0.5 000E 00 -0.I263E-Q6 C,3 P00E Cl -0,15C0E Cl -C.15C0E Cl -0.1 3a3 E-06 

2 K 0.3750F Cl C.9500E 01 0.2500E CO C.ZSOCE CO -0.50COE 00 C.7500E 00 C.7500E 00 0.75C0E 00 -C,15CCE 'ci 0.2250E Cl 


CUMULATIVE STRESS QUANTITIES ♦»»»»<i**»»*»»*»<»***»»*it**»**»»»» * 

ELEMENT EFEECT'iVF ’"’’EFFECTIVE ♦***♦**•*>»** STRESS CENTER ***•**♦*♦**••*•* ***♦***♦*♦***♦, STRESS ******************** 

NG. I.C. CENTER STRESS JX V¥ ll XY XX YY LI 

, a 

tn 

1 20 0_.20J)0E 01 O.AOOOE Q1 0^1333E 01 ~C>6667E 00 -Q.6667E 00 »O.A351E»07 C«»OQOE 01 -0.3863E-C5 C.«2C8E-C5 0«1582E-05 L 

2 lb 0.2000F Cl 0.4000F 01 0.3333E CO C,3333E CO -0.6667E 00 O.IOOOE 01 C.2000E 01 0.2000E 01 -C*39C6E-C5 0,20C0E 01 


I 

IM 


ELEt^M E ~P SUM I NCREMENTAL TCTAL SURFACE **♦* EFFECTIVE PLASTIC STRAINS *** EFFECTIVE CRE EP STRAINS **♦♦ 

NO. I.D. CODE CODE TEMPERATURE TEMPERATURE YIELD SIZE INCREMENTAL SUM INCP. CUMULATIVE INCREMENTAL SUM INCR.' CUMULATIVE “ 

_ 1 _ _ 0 2 C.IOCCE 01 0.6CC0E Cl C.2000E 01 C.IOCOF Cl C.3Q0QF Cl 0.3000E 01 0.0 O.IOCOE Cl C.IOCQE 01 

2 ' 1C ' "0 2 C.IOCOF 01 0.6000E 01 0.2000F 01 O.IOOOE 01 0.3000E 01 0.3C00F 01 0.0 O.IOOOE 01 C.lOflOE fll 




sz-re 


INCPEMEM 


CLKLL A1IVE LCACS ^NC Z-LCACS FCP LCiiC [NCPE>'*ENT 6 

6 


EIEMFK7 I,C. 20 10 

TEMP. C.7000CE Cl C.7CQCGE 01 

7-LOAD C.O C,0 

ITEPA7IVP PPPCP = 0,nC55E-W 


ENC CP LCAC INCPEPEN7 6 


INCPEPEM 6 

HECEAMCAL LOAD CURVE FACTORS = 0.73C0E 01, “o.O 

_CR_EEP TIRE INCREMENT = 0.0 

NO. ELASTIC riEH^TS s' 2, NO, PLASTIC ELEMENTS * 0 

0 ELEME NTS HAVE CHANGED ELASTIC TO PLASTIC, 2 ELEMENTS PLASTIC TO ELASTIC DURING THIS INCREMENT 

SPECIFIEC RAX, NO. STIFFNESS UPDATES * 3, NC, UPDATES PEPFCPMEC * 0 

SPECIPIEC R_AX, NO, ITFPATICNS PER UPDATE = 10, NO, ITEFATICNS PEPFDPMEC SINCE LAST UPDATE = 1 

SPECIFier MAX. UNBALANCEC-FOPCE EFPCP = C.lCCOE-04, ACTUAL EPPCP = 0.17C6E-05 


* ^ **** CU MULA TIV E INTERNAL PDPCES AND C I S PL ACEMENT S ******** 

♦♦ node *• ♦♦******• FORCES •♦•*<r**** •***# DISPLACEMENTS •****•' 

— LcT , __ U V U V 



3C 

-0.1370SC7E-C5 

0.9g«9952F 

GC 

-C.7318404E-06 

0.7299999E 

01 

2 

40 

0.i0860SCE-C5 

0.9999962E 

00 

-0.7711196E-06 

0.7299999E 

01 




-0.20fl2n4E-05 

-0.9999943E 

00 

0.0 

0.0 


4 

20 

0.23669S2E-05 

-0.9999971E 

00 

-0,63fi75C5E-06 

0.0 
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9.1-26 


♦»»»*♦»♦»*»♦»»»»**»:»»»♦»»»»»» ELASTIC STRA INS ***** ***** ** *** ************************ 

ELE^'E^T ThcrM/sL STP4INS *** *4#**i» INCPEKE^TAl ♦♦#**«*i»********* **♦♦**♦****** CUKULfiflVE •»***• •'*♦♦•♦'**••*♦ 


INCPENENTAl 

CLPULATIVE 

XX 

YY 

ZZ 

XY 

XX 


YY 

ZZ 


XY 

0.0 

0.2300E 01 . 

0.0 

0.5982E-C6 

o.c 

-0.2987E-C6 

O.IOOOE 

01 

-0.30COE 

00 -C.300CE 

CO 

0.2153E-C6 

0.0 

0.2300E 01 

0,597aE-06 

-0.A527E-C6 

0.0 

-0i726eE-0T 

C.3500E 

CO 

0.35CCE 

CO -C.30CCE 

CO 

0.6500E 00 


plastic WCPK *♦♦♦ 
INCPEh*ENTAl CUPLLAT 


0.0 _ _ c^95cpe_e.i 

6.0 0.9500E 01 


****♦♦*♦*»»♦»♦♦**♦»**♦♦*«*♦♦**»»»»»* plastic S TP a ins ♦*** **** »»* *** * »» *»* * *»* *?*•***••_ 

*,^m*****ni^*^* IMCPEPENTAt ^ ••♦♦♦**#**>•*' cuf'ULATIVE ' *»•***'»•*****•**•• 

XX 5fY U XV XX YY ZZ XV 


C.3000E 01 -0.15COE 01 -C.15CCE Cl -0.1383E-06 


C.7500E CO 0.75C0E CO -C.15CCP Cl 0*2250E Ol 



_ _ *♦♦♦**»»*»♦*»♦*♦**«»♦*»*■*♦*»» CU »*UL A T I V E ST P ES S CUANT I TIPS *» »»* *** ***» *»»**♦**»******»»*♦♦ * 

ELEPE^^ EFFECTIVE ^EFFECTIVE ************ STPESS CENTEP **************** ■■#*******♦***•** STRESS ****♦*♦**♦*•****•••• 

NC. r.C. CENTE R STPESS XJ< YV Lt XY XX YY ZZ XV 

1 20 Oa^gOOE 01 0.2000E 01 0.1333E Cl -0.6667E 00 -Q.6667E QC -0.4351E-.C7 C.2QOOE Cl 0.A4C3E-C5 0.6898F -05 ^3^2E-06 

2 10 o'.200bE bl 0.2000E 01 0.3333E 00 C.3333E 00 -0.6667E OC C.IOCOE 01 C. lODOE 01 C. IDCOE 01 -C,ie32E-C5 C.IOCOE 01 
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9.1-27 

I 


CUHLLATIVE T HEP>*OL LCflCS ANC Z-LOADS f^CP t O AC INCfiEKCN T J 

INCPEMFKT 7 


EIEMeKt !.C. ’20 10 

TEMP. C.54CC0f C1_C.5400J)E 01. 

Z-LCftD C.C C.O 


I TEPATI vr EFPCP « 


tTEPAT IV E gppnp = 


ITERATIVE EPPCP - 


ITERATIVE EPRCF - 
ITERATIVE EPRCR * 
t T EOAT l VE EPPCP « 
ITERATIVE FPPCR » 
ITERATIVE EPRtjP e 
iterative ERPCP = 


|TE_PATI VL-EMCf = 

rTEP“ATivE ERPCP - 


0.42622E 00 
A 3J171P 00 


0.23360E OC 
fl aiC 15J 
0.12ACeF CO 


CL..S.CiSlfrLQl 

0.<E250E-01 

Q.4eE30E~01 

0.355S7E-01 


0.26137L-0J. 

Q.12619E-C2 


ITE RATIV E E PPCP <» 0. 11159E’0 3 
ITERATIVE ERPCP » 0.12A30E-C5 
ITERATIVE ERPCP a 0.2431SE-C6 


! 


I 


ENC CF LCAC 1^CPERE^T 7 


INCREMENT t" ''' ' 

^ECtiAMCAL LCAC CUP VE FACTORS « 0.60C0 E 01 « C.C 

CREEP TIME INCREMENT * 0.0 ' ^ 

NC, ELASTIC ELEMENTS = C. NQ. PLASTJC ELEMENTS ^ 2 

2 ELEMENTS HAVE CHANGED ELASTIC TO PLASTIC, 0 ELEMENTS PLASTIC TC ELASTIC DUPING THIS INCREMENT 

MAX,_NO, STIP FNESS UPDATES « 3, NC. UPDATES PEPFCPMEC g 1 

SPECIFieC MAX. NO- ITERATIONS PER UPDATE « 10, NG. ITERATIONS PEPFCPMEC SINCE LAST iJPCATB = 4 

S P E C I F I E C MAX, UNBALANC ED~FCPCE ERPCP - O.LOCOE-04, ACTUAL ERROR = 0 . 3431E-06 



******* CUMULATI VE INTERNAL FORCES AND D I SPL ACF MENTS ******** 

*♦ NODE ** *♦♦♦***** forces ****♦♦*** ***** DISPLACEMENTS ****** 

. NC, I U V U, y 

1 30 >Q.41723 2EF-06 -0. 35 16E74E-05 “0.6S71468E-06 C.6C000C0E 01 

2 40 6.49875e5E-06 -0.2861023E-05 -0, S3 36210E- C5 C.6COOOOOEOI 

. 3 1C ,.-0.l36lCC4E^.6__0,405 3J16E>C5 0.0 0.0 

4 20 0.54576S6E-07 0. 23245EIF- C5 -0.698B167E-05 0.6 


8,2-1 *6 


»4»»^-»»»»**»»»*****»*»»***»»<*» »*»»* glASTlC STRAINS ****♦«»♦♦*♦♦*♦» 

eLEMFM *** THERV4L STRAINS ♦♦* ♦**♦♦»*♦•**** ' INCREMENTAL ***************** ************* ' CUMUtAt’lVE "******************' 


NC.. l.C. incremental 

cumulative . 

XX 


YY 


zz 

XY 


XX YY 

ll 

XV 

1 20 -0.3CCCE 

cc 

0.2003E 

Oi 

-O.IOOOE 


0.3000E 

00 

0. 300CE 

00 -0.8922E- 

06 

0.4172E-C5 -0.78CaE-05 

0.13UE-C5 

-0.6769E-O6 

2 10 -0.3000E 

cc 

0.20C0E 

01 

-C.3500E 

CO 

-C.3500E 

00 

0.3000E 

00 -0.65CCE 

00 

0. l7BaE-06 -b.476aE-06 

O.U92E-06 

0,4530E“b5 


»»**»»»»*»»»»**»»»)»♦*♦»»»»♦»»»»»*»» * PLASTIC STRAINS *************************************** 

element *♦*** PLASTIC hCRK **** ************* INCREMENTAL ***************** ' ****♦**♦•*♦•♦ CUMULATIVE ♦***♦*♦•*•♦*♦***♦* 

NC» r.C. incr emental cumulative XX VY U XV XX VY, ZJ XY 

1 2C 0.0 _ _0.9500E 31 - 0.2264E-05 0.U3ZE-05 0.1132E-05 -0.B163E-12 C.3Q0QE CL -Q .LSOQE 01 --<S «I5QCE 01 -0«1383 E-06 

2 IC 0.0 ' C^SSOOEOl -C.3509E-C6 -C.3509E-06 0.7C19E-06 -0.1053E-C5 0.7500E 00 0.7500E 00 -0.15C0E Cl 0.2250E 01 


_ _ ***************************** CUMUL ATI VE ST PES5 CUANTITIES ********************************* 

element”’ effective effective ************ STRESS CENTER ♦♦♦*♦♦♦♦**♦•♦**♦ *************** STRESS ******************** 

KC. l .C. CENTER STRESS iiS H U 51 XX YY U XY 

1 2_p 0.2000E 01 Q. 7fl67E-05 Q.1333E 01 -0.6667E 0 0 -0.6667E 00 -0.43S1E-C7 “C. 333BE-0S 0.6298E~06 O.56b5E-0S -0.5207E-06 

2 IC 'c.2066f ci C.A907E-05 C.3333E CO C.3333E C0-O.6667E 00 O.IOOOE 01 -0.29AAE-05 -0. 1729E-05 -C.33576-C6 -0.2623E-C5 


ELEMENT E-P SUM INCRIMENIAL TCTAL SURFACE **♦* EFFECTIVE PLASTIC STPAINS *** ***** EFFE CTIVE CREEP STRAINS **** 

NO. I.r. CODE'CCDE tERPERATURE TEMPfRATUFE YIELO SIZE INCREMENTAL SUM INCP. CUMULATIVE INCREMENTAL SU’m INCP. CUMULATIVE 
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8 


CL MULA Tt VE TH E P(*AL LG i >CS JNC Z-LCOCS FCP LC^t INCFEKEM 

INCPEWEM 8 

"TleVfTt' iTcT 20 To 

TEMP. C.8Q000F Cl C.80000E 01 

Z-LC8C -C.60000E CC-C. 60000? 00 

ITEP»TIVE PPPCP i ”0 Ts5£20E 00 

ITERAT IVE EPPCP - Q.76S52E CO 

ITEPATUE EPPCP * 0.61256E 00 , 1 

ITERATIVE EPROP = 0.26E32E_00 

ITERATIVE EPPCP = 0.10396E 00 

_ ITERATIVE .ERROR - 0.29551E- C1 

ITERATIVE ERROR » O.S IA52E-C2 i 

ITERATIV E g PPCP « 0.275380-02 ... 

iterative error = 0,e36C4E-03 

ITERATIVE ER_RCP = 0.2 S2a9E~03 ^ 

ITERATIVE ERROR = 6TlCe76E-04 

ITERATIVE ERROR f_ 0. i295 9E-C5 ^ 

ft ER AT ivf ERROR * 0.59 506E-06 



INCPE_HEKT_ _ 8__ 

PECEAKICAL lOAc 'CUR VE ' FACTORS = 0.4CCOE 01, C.O 

CREEP tire INCREMENT * 0.0 _ ^ 

NC. ELAftlC ELEMENTS » ' 'O, NO. PLASTIC ELEMENTS = 2 

0 ELEMENTS HAVE CHANGED ELASTIC TO PLASTIC, 0 ELEMENTS PLASTIC TO ELASTIC DURING THIS INCREMENT 

SPECIEIFC MAX, NO. STIFFNESS UPCATFS = 3, NC. UPDATES PERFCPMEC * I 

SPFCIFIEC MAX_. NO. ITERATIONS P EP U P CATE = 10, NO. ITERATIONS PERFCPMEC SINCE L AST U PD ATE = 

SPECIFIED MAX.' UN BALANCE D-FCPCE ERRCR = 0.10 COE- 04, ACTUAL ERROR » 0.59516-06 



♦♦ NODE ** 
- N L> L,X^. 

1 3 0_ 

2 40 
_ 3_ - 1C 

4 20 


<,♦»» »«» CUMULATIVE INTERNAL 
«*«»«**»* FORCES *♦*♦*♦*** 


=CPCES AND DISPLACEMENTS *******w 
***** DISPLACEMENTS ♦♦♦•** 


C.2980233F-06 -Q. 3L2S 0 35E 00 
0.75961916-07 -0.3125043E 00 
0.37 4C4 5 IE-06 Q.3125044E 00 
0.1519833E-06 0.3125035F OC 


-0.1033e93E-C5 O.4CO O0O OE 0 1 
-C.600CC21F 00 C,4C06 cC06 01 

0.0 OiQ 

-C.6000C34E 00 0.0 


2-992iL- 





******* 4 *4 **««%«»«*» «!»****«< ELASTIC STRAINS ********** «*4*********»***«********«*4* 
ELEMEKT *** THEPJ*AL STRAINS *•* ************* ' fNCPE KENT AL ***************** ************* CUI'ULATIVE ****4************* 



I.C. INC REM ENTAl 

CUMULATIVE 

XX 


YY 


ll 


XY 


XX 


YY 


It 

XY 

1 

2C -C.ICCCE Cl 

0. lOOOE 

01 

-0.5000E 

00 

Q.1500E 

00 

0.1500E 

00 

0.2709E- 

06 

-G.50Q0E 

00 

0.1500E 

00 

0.1500E 

00 -0.<i059E-06 

2 

10 -o', lOOOF 01 

o.idooE 

01 

-0.1750E 

CO 

-0.175CE 

CO 

0.I500E 

00 

-0*3250E 

00 

-0.1750E 

00 

-0.175CE 

00 

C.15CCE 

CO -0.3250E CO 


****«»**4****** * *4**»»»******»»»*»*« pl astic STRAINS ********** ******** ********** **•**♦***•• 

ELEI^ENT *****"’>LASf 1C WCPK ***♦ ************* INCPEf'FNTAL ***************** ************* cTJi^LATI VE ******* •******««**'" 

KC. !.0« INCREMENTAL CUMULATIVE XX YY XY XX YY ll XY 

. I 2C 0.1563F 00 Q. 9656E 01 ~Q«5Q0 0E 00 a»2S0QE 00 O,25Q0E QC -0.92UE- 07 C - 2500E 01 -Q. 1250E 01 ~0.1 2 50E 01 “0.23Q4E-Q6 

2 10 0,1563P CC C.965&E 01 -C.125CE C0-0.12S0E CO 0.2500E 00 -0.3750E 00 0.6250E 00 0.6250E 00 -C.125CE Cl o'-iaT^E Ol 


»»4**»*»******4*»»* 4 **♦»»*♦** COPULATIVE STRESS CiMNTITIES •***♦•***•** ••**44 •*♦*•****♦*♦*•• 

ELEMENT EFFECTIVE EFFECTIVE *****4****** STRESS CENTER *****4**»******4 *4*******4***** STRESS ************* «'*^*«* 

NC. I .C. C ENTER STRESS XX YY ZZ XY XX YY IZ XV 

1 20 0.,150^E Ci 0» 6250E 00 0«10Q0E C l - 0.50Q0E 00 ~0«50 00E 00 -0.1049E- 06 -0.6250E 00 0,a637E-07 C.22C1E'05 -0. 3903 E- 06 

2 10 0.1500E 01 C.6250E 00 0.2500E CO C.25CCE CO -0.5000E 00 0.7500E 00 -*0,3125E 00 -b.3125E 00 -0.3228E-05 -0-3125E 00 



eiE ljEM E-P SUM INCRERENTAL TOTAL SURFACE ** »* EFFECTIVE PLASTIC S TRAIN S *»* * 4 * 4 * EFFECTIVE CREEP STRAINS 

NO. I.I3, CODE CODE TERPERATURE TEMPERATURE YIELO SIZE INCREMENTAL SUM INCR. CUMULATIVE INCREMENTAL SUM I NCR. CUMULATIVE” 
1 Q 2 C.26CCE 01 0. SCOPE 01 Q.2125E 01 C.500CE CO 0.3500E 01 0.2500E 01 0-0 C.iOCOE Cl G-IOOOE 01 
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lE-l-6 


9 


INCFfMEM 9 

ELErETT”T7c~ ' 2^ 10 

TEMP* 0.90000E 01 C. 900008 0 1 

Z-LCAC *C.20b00E CC-C.2COCOf OC 

ITEPATIVE EP^P' = “o.e^^iSE 00 
ITEPAT IV E ERRCP g 0.74A55E CG 
ITfPftTIVE EPPCP s 0.4G913E 00 

. tJePATjVE .ee»CP 0. I4128g-Ql 

ITERATIVE EPPCP « 0,tS163E-02 

.1 TFPATl^ E PRPP g 0.12C6BE-02 _ 

ITEPATIVE EPRCP s 0,29t64E-03 
ITEPATIVE EPRCC = 0.t41B5E‘04 

ITERATIVE 6PBGR = 0.14244E-04 

Ilf PATtv e EPBCP = 0.41961 E -05 

ITEPATIVE ERRCR « 0.fcn26E«05 


END CF tCAC IKCPEMFNT 9 


_JNCF£MEAT 9 

PECHAMCAl LOAC CURVE FACTORS = 0.1500E 01, oTo ^ ' 

-CREEP T I Mf _ 1 N CPEM ENT » _ O.I QOQE 02 

NO. ELASTIC eIemeMS = 0, NO. PLASTIC ELEMENTS = 2 

— 0 - - f -C.F P E NI5 -HAVE .C HANGED ELASTIC TC PlASTlCt 0 ELEMENTS P LASTIC TO E LASTIC DUPING THIS INCREMENT 

SPECIFIED MAX. NO, STIFFNESS UPDATES * 3, NC. UPDATES PEPFCRPEC * 1 ' 

—SPECIFIED FAX, NC. ITFPATICNS PER UPDATE = 10, NQ« ITERATIONS PERFORMED SINCE LAST UPDATE = i 

SPECIFIED MAX. UNPALANCED-FDPCE ERROR = O.lOOOE-04, ACTUAL ERRCR > 0,61736-05 " 


♦**♦**♦ CUMULATIVE INTERNAL FCRCES AND DISPLACEMENTS »*»♦*♦*» 

NODE *• ********* FORCES ♦**♦**•** ***** DISPLACEMENTS •»*•** 


NO. 1 

:.c. 

L 

V 


u 

V 


X 

3C 

-0.2205371E-C5 

-0.6250C15E 

cc 

-C.3952121E-C5 

C. 15Q0CC0E 

01 

i 

40 

0,14398696-05 

-0,6250042F 

00 

-0.2000027F 00 

C, 1500066e 

01 

.3 

IC 

0..25761C2E-C5 

0.625CC3CE 

cc 

C.O 

0.0 


4 

20 

-0.l8i06C06-C5 

0.6250029E 

CO 

-C.200005CE 00 

C.O 
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**»*<»»** gLA STI C STP/iI_N^ it^jm ^^m********^*** * * **** ***** * ****** *** 

ELEMEKT ♦**” THERMAL STRAINS *** ************* INCREMENTAL »•»*#••********♦*■ CLMULATIVE" ****************** 


NC. 1 

I.C. INCREMENTAL 

CUMULATIVE 

XX 


VY 


n 


XY 

XX 


YY 


ii 

XV 

1 

20 -0.5000E OC 

O.SQOOe 00 . 

-C.5CC0E 

00 

0.1500E 

00 

0. 1500E 

00 

0.1456E-C? 

-C.IOCOE 

01 

0.3000E 

00 

C.300GE 

00 -0.3914E-06 

.. 2 

10 -O.SOOO'f 00 

0,50006 00 

-0.1750E 

00 

-C. 1750F 

00 

0.150CE 

CO 

-Ci3250E CC 

-O. 356 OE 

00 

-O.35C0E 

do 

0.30CCE 

CC -0,65 COE do 


************************************ plaST 1 C STRAINS *************************************** 

ELEMENT PLASTIC WORN **** *♦***•♦♦**♦♦* INCREMENTAL *♦•*•*•**♦**♦***♦ ************* CUMULATIVE **♦#**••*♦•*♦*♦*•♦ 


NC. 

I.C. 

INCREMENTAL 

CUMULATIVE 

XX 

YY 


IL 

XY 

XX 


VY 

ll 


XV 

1 

20 

0.4686E OC 

0. 1012E 02 

-C.50C0E 00 

0.2500E 

CO 

0.2500E 

00 -0.5287E-07 

C.2000E 

01 

-0. lOCOE 

Cl -c.iococ 

Ql 

-0,28336-06 

2 

10 

0.A688E dc 

0.1012E 02 

-C.1250E CO 

-0.1250E 

00 

C.250CE 

00 -Q.3750E CO 

C.5000E 

00 

0. 50 COE 

do -c.iocct 

Ql 

O.L'SQOE Ql 


»*♦»»»»♦*»♦»**<>♦»♦»»*♦»*»♦**»»»*»»»* CR EEP S TRAINS **♦*♦*♦*****♦♦*♦**»***♦♦*♦*****♦♦•***♦* 

ELEMENT ♦♦***’* CREEP WORK ***** •***•♦••#***♦ INCREMENTAL «***»*«««*«** CUMULATIVE *•*•*••*•♦♦♦*♦**•* 

liL* l.-C. INCREMENTAL _ CU Mm A T_LV E U il U Xt XX Tt It XV 


1_ 20 0.9375E 00 0.3 6 88£ 01 -C.ICOO E C i C.500CE CO 0 iSO 0^ 00 -0.4683E-06 0.238AE-06 0.2 3 6^E-Q« 

2 1C 0.9375E CC 0.3688E 01 -C.2500IE CO -C.2500E CO 0.50C0E CO “C.750CE 00 C.4768E-Q6 O.1O73E-0S 


-0.A 768E-06 -0,A791E~06 T 
-0.1550E-05 -0.5960E-07 ^ 

(S5 

o> 

(Ti — 

I 

IS? 


»»»»»»»»♦»♦»»»»♦»»»»»»»»»»»»» CUMULATIVE STRES S QUANT ITIES ♦»»»»»»»»»»»»» »»**»»*»*♦»♦»»»♦»»» 



element 

EFFECTIVE 

EFFECTIVE 

************ STRESS 

CENTER **************** 

*************** STRESS 

*«*****«« *41 


NC. 

1 .C. 

CENTER 


STRESS 


XX 

YY 

ll 


XY 

XX 

YY 

ll XY 


1 20 

O.IOOOE 

01 

C.1250E 

01 

0.6667E 

CO -0.3333E 

GO -0.3333E 

00 

-0.X402E-O6 

-0.1250E 

01 -0.3245E-C5 

C.1969E-C6 -C.3763E-06 


2 

10 

C.IOOOE 

01 

0.1250E 

01 

0.1667E 

00 0.1667E 

00 -0.3333E 

CC 

G.500CE 00 

-C.6250E 

00 -0.5250E 00 - 

0.1036E-05 -0.6250E 00 


_ ELtMENl E*P SUM IKCPE^ENT AL TCT4L SUPfACE » **» EFFECTI VE P LASTIC STRAINS *♦* ***** EFFECTIVE CFEEF STRAINS 

NO, UC. CQCE CQOE TEFPERATUFE TEMPERATURE YIELD SUE INCREMENTAL SUM INCP, CUMULATIVE INCREMENTAL SUM INCP* CtMLLATIV 6 ~ 

1 20 0 2, C,1000E 01 0«90COE 01 0,2250E Cl 0,5CC£E CO C^4QC0E Cl Q>2CQQ£ Ql C*10a0E 01 C*2000E G1 0,7304E-Q6 

2 10 0 2 C,IOCOF 01 C,90COE 01 0,225CE 01 O.SOOOf GO O.AOOOE 01 0,2000E Cl O.IOOOF Cl C.20CCE Cl 0*15fi9E-05 




CIMCLATIVE THE P^Al LCfltS j>NE Z-LCACS FCP LCAC IKCPEPEM 10 

rNCPEME^7 10 

eTeV eSt ■ ] .n , 20' To " ^ ^ 

levp. O.IOOOOE 02 C.IOOOOE ^2 . 

Z-LCAO C.C C.6 ’ ■■ ■• . 


ITERATIVE 

-tl£PATLVE_ 


EPPCR 
PPPCF_ 
ITERATIVE EPPCR 
ITERATIVE £P«CP 
iterative' eprcp 

ITJ.PATIVF EPKE 
ITERATIVE EPPCR 


C.936C6E CO 
0«S320ZE 00 


0.E4505E 00 
0,i0C3_4E 00 
0, 61599 E~Ol 
0 .56229 E-C1 
O,A5eiOE-0i 


ITERATIVE ERPCP = 0.25E34S-01 

. . ^lJEPAJli.L.EPRj:P * 0.1912AE-01 

ITERATIVE EPPCR = 0. 

ITERATIVE EPRCP » 0.H579E-06 


_ 



CJ 

cn 

~Z- -- j. 

N C C F 

LCAC INCREMFNT 10 

fV3 




tr> 

t 

ro 

INCREMENT 10 

MECiiA.NICAL LQAC CURVE FACTCPS = 

0.365CE 

Cl. 0.0 


CREEP TIME INCREMENT * 0.0 

_N0. ELASTIC ELEMENTS = _2j. NO. 

PLASTIC 

ELEMENTS = C 


U fcLhP'fcNIS HAVE CHANGED ELASTIC TO PLASTIC* 2 ELEMENTS PLASTIC TC ELASTIC DURING THIS INCREMENT 

SPECIR.IEC MAX,_ NO.. STIFFNESS UPCATES = 3, NC. UPDATES PERFCPMEC « 1 


SPEcr^iec MAX. no, itepaticns pep update 

SPECIPIEC max, UNEALANCEC-FCRCE EPRCP = 

= 10, NC, ITERATIONS PERFCPMEC SINCE LAST UPDATE =2 

O.lOGOE-04. ACTUAL ERROR * C.U58E-06 













*»..*!*** CUMULATIVE INTERNAL 

FCflCES AND CISPLACEMENTS 


♦* NODE 

NC. 1 

* + 

. c 

i»»*i^**+»* FORCES **♦*♦**•♦ 
U V 

***** DISPLACEMENTS ****** 
U V 


1 

30 

0,17881396-06 0.A9S9S76E 00 

-0.S417178E-C6 

0.3650000E 01 


2 

3 

AO 

1C 

-0.1528397E-C6 O.<t999906E 00 

C.IC812C7E-07 -0.A99997CF 00 

0.^tl72325F-06 

0.0 

0.3650CCOE 01 
0.0 


<1 

20 

-0,36786336-07 -C.SOOCCCCE CO 

-C.U92093E-C5 

0.0 ■ 

• - 




9.1-34 



»»*»»»»» »»»!>»♦»»»♦»»»»♦♦#■ «*»»»**♦*»* EL /> STIC STP M NS _ ** ** **** * ? ? * ** * * ** ***** 

ELEHEM •** thermal STRAINS *** ************* INCREMENTAL ***************** **♦*♦*%♦****• CUMULATIVE •**♦♦**»♦*****•♦*• 




1 2 C 0 . 65 CCF CC Q.U 5 QE 01 . 0 . 15 QQE Q 1 -O.A 5 JOE 00 - 0 « 4500 E 00 Q, 3576 E-C 6 C. 5000 E 00 - 0 . 150 0 E C O -C.I 5 COE CO - 0 . 33 T 6 E -07 

2 10 0 . 6500 E GC 0 .U 50 E 01 C. 5250 E CO 0 . 5250 F CO - 0 . 4500 E 00 C* 9750 E 00 0 . 1750 E 00 0 . 1750 P 00 -C-I 5 CCE 00 d. 3250 E CO 


I 

t 


! 


»»♦********♦♦»♦»»♦»*♦*»*»**♦****»*•* PLASTIC STRAINS ♦**♦***♦*♦*♦** ****** ***** *»* • M * *♦**♦*• 

ELEMENT ***** PLASTIC MC'PK **** ************* INCREMENTAL ***************** ************* CUMULATIVE ****************** 




INCFEMENTAL 

CUMULATIVE 

XX 

YY 

zz 

XY 

XX 


YY 

ll 


XV 


1 

2 C 

0 . c 

0 . 1012 E 02 

0.0 

0.0 

0.0 

O.Q 

Q.ZOOOE 

Cl 

- 0 . lOCOE 

01 -C.IOOCE 

01 

- 0 . 2833 E -06 

i 

1 

2 

10 

0.0 

0 . 1012 E 02 

C.O 

C.O 

0.0 

o.d 

C. 5000 E 

00 

O. 5 OC 0 E 

CO -C.ICCCE 

ci 

C.lSCOE Cl 

1 


! 


EFFECTIVE EFFECTIVE 
CENTER STRESS 


* *»**»*»******»♦»***»»♦**♦*** CUMULATIVE 5 TPES S CUANTI TI E^ ***** *^* * * »♦* O • ^ ^ *♦*♦* 

************ stress CENTER **************** ♦♦***♦•**•*♦**♦ STRESS *******************o' 

XX o U o XX VY It XV 


L 2 .Q_P.IOLPOE 01 0 , 10 QQE 01 0 « 66 67 E 00 ~Q» 333 3 F 00 “ 0 « 3333 E ^ j^,_lAq 2 E:iC 6 _ lOOOE _Ql _- 0 i, 2 U 3 E -01 0 .<> 12 CE-C 5 - 0 . 5 I 95 E ~07 

2 "id o.i’odoi 01 O.IOOOF Ol C, 1667 E CO C. 1667 E 00 - 0 . 3333 E OC 0 . 5000 E 00 O.SOOOE 00 O.SOOOE 00 - 0 . 25 A 0 E-C 5 0 , 50 CO£ 00 


ELEMENT E-P SUM INCREMENT AL TCTAl SURFACE **** EFFECTIVE PLASTIC S TRAINS *** ***** EFFECTIVE CRE EP STRA IN S *♦* » 

NO, 1 , 0 . COnE CCOE TEMPERATURE TEMPERATURE YIELC SUE INCREMENTAL SUM INCP. CUMULATIVE INCREMENTAL 'SUM INCP. CUMULATIVE 

1 20 - 1 -2 C.IOCCF 01 C.IOCOE 02 0 . 2250 E Cl Q.O C. 4000 E 01 0 . 2 C 00 E 01 0.0 C. 2 CCCE Cl C. 7304 E -06 

2 1 C -I -2 C.IOCOE 01 O.IOOOE 02 0 . 225 CF Cl 0.0 C. 4000 E 01 0 , 2000 ? 01 0.0 C, 20 C 0 f Cl 0 . 1589 E -05 


05-17266-2 



SE-L‘6 


CL >*ULATI Ve THEPyAL LCACS ONC Z-LCACS FC P LGAC INCPEHE NT IJ^ 

INCPEMEM 11 

£LEPfTt 20 To 

TEMP. C.45000F 01 C.^SOCO E 01 . 

Z-ir/iO -0*30000E CC-C.300COE OC 

ITEPATIVE PPRCP = ' ’orilASTTTo 

HERAT I VE EPRCP * 0.S3666E 00 

ITEp'at'i've EPPGP * 0.5«S<|9E 00 

ITEPATIVE EPRCP * 0.1 5193E CO 

ITEPATIVE EPPCR = 0.1C629E CO 

ITERATIVE EPRCP * 0,{2 319E-0 1 

ITEPATIVE EPP’CP » O.^C663E-Ol “ 

TTEPAT I VE EPPCP - 0,24767E-01 

ITEPAT'tVE EPRCP * 0.15738E-01 

JTEPATIVE EPPCP » O.^J163E*i02 

ITEPATIVF EPRCP = '0.17456E-06 ' ” 

ITERATIVE EPRCP = 0.2S0A IE-C5 



INCPEPEKT ll 

MECHANICAL LQAC CURVE FACTORS * 0.5500E 01, 0,0 

CREEP tlME'TNCPEHFNT = O.IOOOE 02 

ilO. ELA^TC ELEMENTS = 0, NO. PLASTIC ELEMENTS = 2 

"2 'ELEMENfS HAVE CHANGED ELASTIC TO PLASTIC, 0 ELEMENTS PLASTIC TC ELASTIC DUPING THIS INCREMENT 

SPECIFIEC MAX. NO, STIFFNESS UPDATES = 3, NC, UPDATES PEPFOPMEC = I 

SPECIFIED MAX. NO. ITERATIONS PEP UPDATE * 10, NO. ITEPATICNS PERFOPMEC SINCE LAST UPDATE = 2 

SPECIFIED MAX. UNBALANCED-FCRCE EPPCP = O.lCCOE-04, ACTUAL ERROR = 0.250AE-OS 



1 3C -C.22Q5371F-05 0.162<i9q9E 01 -C. 1 A 769 1 2E- C5 C .550 CCC0E 01 

2 AO 0.19148E8E-C5 0.l62<t997E 01 -C,2999905F CO C.5500CCOF 01 

3 , . . 10. 0., 4 1 3 8 4 2,3 E -0 5, , - 0 , 1 6 g A 99fiE 01 0*0 0.0 

A 2C -C.3e4792CE-Q5 -0.1629995F 01 -C.3C00008F 00 0.0 


D5-1 7266-2 





D5-U266-2 



aAiitfin/vni *dD'ji .^ns n?i'jdc*3aDNi SAUviTwriD 'aoNi wr>s iffiNcj^aaoNi 37is aiaiA 3ar>iva3di><3i 3afuvaad(^3i 3033 3ao3 *3*i 
• •♦♦ Slstivair d33(J3 3AI1333d3 '*** SVl'9tilS 3ilSVld'3AU)33'33 •*** 3D»ddr>S IVlOi IViNlHWadDNI WHS d-3 1»43W313 



10 35Z91*0 

S3-3US1*3 

- 13 35Z91*0 

13 

35291*0 

00 30005*0 

DO 

at£f£*0- 00 

37991* 

0 00 37991*0 

to 30SZE*D 

TO 30001*0 

01 

Z 

90-39ZZ.£*0- 

50-35971*3 

53-3eZ£7*0- 

ID 

30SZE*3 

9o-azo^i*o- 

30 

3££€£*0- 00 

ieeee* 

0- 00 37999*0 

10 305Z£*0 

to 30001*0 

02 

T 

AX 

21 

AA 


XX 

AX 


U 

AA 

XX 

SS3diS 

d31N33 

•3*1 


t *« t* tl t «t 1*111 ttttt 


11111111*1111111 d 3 IN 3 3 

SS3diS 

111111*11111 

3AI133333 

3AU33333 

iX3ri313 


«*«i^*«**«*«K* *************** **««* S3 1 II iJ^vnb ssBdiS 3 a i iv *»**»***********♦♦*♦**<(♦**•*♦ 


00 30Dfii*0 00 33005*3- 00 300Se*0 03 300SE*0 00 30051*0 DO 30005*0- CO 3D05Z*0 00 30052*0 tO 3£T8S*0 10 JSZIZ'O 01 Z 

90-3«?lZ/.*0- 03’ 30305*3- bo 30005*0-13 30001*3 90-3eZ^Z*b- 00 30005*0- 00 30005*0- ID 30001*3 10 3C l®5 *0’ 'l 0 35 Zl Z *0 ' OZ I 

AX n n 7% ix n n th 3^11^10^03 mNii« aoNr '•a*i *w 

«••******«««« Hit 41 3Al iyir>d03 *T***** <i»»»*»»»»« 1XlN3 /<3d3Nl ♦»»»»»»»♦»♦»* >taOW d33a0 1X3^(313 

«***#^«*#** 1* #«*««« **.*««« SNlVHiS 33383 miiniiiiiiiinmmiitimimiii 



10 30051*0 13 30301*3- 00 30305*0 00 30005*3 0*0 0*0 0*0 0*0 ZO 3Z101 *0 0*0 01 Z 

90-3E€BZ*0- 13 30301*3- 13 303’01*0- 10 30002*3 0*0 0*0 0*0 0*3 ZO 3Z101*0 ' 0*0' OZ I 


AX“’’ 77 AA Xi AX n AA XX 3AI iVindID 1Vi,X3H333NI *3*1 

3Mi7 ind03 1VlN3iX3d3Nl *»«★ Xd OH OJ^SVld n*** iX3H313 

iiiiiiiiiliilliiiiiiliiimiiiii SNWdiS OliSVId 


00 30359*0 

33 

3D30£*3- 00 

aoosE’o 

00 

3005£*3 

00 305Z£*0 

00 

30051*0- 30 

30S7TO 

CO 

30S7TD 

TO d03ST*O 

03 JDOSE’O 

31 

?. 

90-31591*0- 

33 

3b00£’*’3"- ’Ob 

3030C*0- 

10 

30001*3’ 

90-aE’sii*b- 

00 

30051*0- 00 

3D05TD- 

03 

33005*0 

■ TO 30051 *0 

’0 0 300S£*0 

OZ 

1 

AX 


u 

AA 


XX 

AX 


u 

Kk 


XX 

3AIiVl1dfl3 

1VlN3W3d3NI 

’•r“l ■ 

■ *0'N" 


11 It «« «i « 3AIAVin^l^^ ii tittitiitti i» tiitii iiiiiiiii 1tfiN3d3 33Xl iimitiiti ii *il SNiJdiS 1?aa3Hi ii« iX3W313 

♦ ittiiti* iiitlititititltitititiitii’tii* SfxiffaiS 3lISV13 iiiititiiiiiiitiiiiiiiiiitiitttitiii 


9.1-36 



4e-i*6 


C t .I? tLii T IV E THEPVaL LCACS Z-LPADS P CR L 0 AC_I NCP Ef^ENT 


12 


incpevem 


12 


ELEf'EK'T I.C. 20 10 

TEf^P. C.,12C00E G2 C,12000E_02 • 

2-L040 O.C C.O 


ITEPATIve 
ITE PATIVE 
I TE PAT IVE 
I TEPATI VE 
riEPATIVE 
I TEPATIVE 
ITEPATI VE 
I TEPAT I V E 
ITERATIVE 


ERRCP 
EilP.fP 
6 PROP 
EPRCR 
ERP'CP 
5PRCP 
EPRCP 
_EPPCP 
EPPC’P 


0.S3777E 00 
Q «ece30E 00 
0.67C5<*E CO 
0. 59i eO F~0 1 
C.E3E50E-C2 
0,ES916 E-03 
6.6S355E-04 
0. 76S21E-05 
0.5C321E-G5 


END CE LCAC INCPBPENT 12 


INCREMENT 12 

'MECPASIcU 'LOA'c CUWe FACTCRS » 0.65 50E Cl, G.O " 

CRE EP TIPE IN CREMENT » 0,0 

NG. ELASTIC ELEMENTS = Ct NO. PLASTIC ELEMENTS = 2 ^ 

0 ELEMENTS HAVE CHA NGEC ELASTIC TO PLASTI C, 0 ELEMENTS PLASTIC TC ELASTIC OLRING THIS INCREMENT 

SPECIFIEC MAX. NO. STIFFNESS UPDATES * 3, NC. UPDATES PEPFCRMEC = ' 0 

SPECIFIEC MAX. NO. ITEPATICNS PER UPDATE = 10, NO. ITERATIONS PEPFCRMEC SINCE. LAST UPDATE = 9 

SPECIFIED MAX. UNEaIanCED-FCPCE EPRCP » O.lOOOE-04, ACTUAL ERROR * 0.5032E-05 ' ' 


. ••.***** CU MULATIVE INTERN A L FC PC ES A ND C 1 SPL AC EMENT S **♦**♦*♦ 

•• NODE ♦* i******** FORCES *•***♦♦** »♦*•* DISPLACEMENTS •♦•*** 


NC. 

l,p. 

y 

y 


u 

V 





. .1 

. 30 

. “C,_1.72.853AE-Q1. 

0. 1937AS5E 

Cl 

-C.A1725G7E-C6 

0,6549999E 

01 




2 

AO 

0.21625T5E-C5 

0.1937A95E 

01 

C.125169eE-C5 

C,6549999F 

01 




... 3 

10 

-C,g9C69ECE-C6 

-0. 1937A95E 

.01 

0.0 

0.0 





4 

20 

0.5566571E-C6 

-0. 1937A96E 

Cl 

-0.65565UE-C6 

O.C 






05-17266-2 



C'*D- 


TO 300CTO 13 ^000£*0 


TO 30001*3 13 3030C*0 

3A11P13413 *d3VI ^ns 


0*0 


*«** SNItfdlS d33dO 3MJl33d33 «*«** 


TO 300SZ*r TO 300S»*0 33 33D0S*0 
TO 300SZ*0 TO 300S«^*0 00 30005*0 
_^Ji.ynwn3 *if3NI WfTS Ttri?^3W3d3Nl 
*•* SNIVSiS DliSVld 3M133d33 **** 


ID 351£2*0 
TO 3StEZ*0 
37IS 0T3IA 
rSVdilOS 


ZO 30071*0 TO 30DSi*0 2 0 01 Z 

ZO 300ZTO TO 30051*0 Z 0 0? ! 

3(mjyd3dw(3A daniffd acIWdi 30J3 3J0D *3‘1 *on 
IV iOl WTNTd3d3^^I tins d-3 iN3H3T3 


CM 


TO 3i.£6T*0 60-»i790»*0 TO 38 E&t*0 TO 30€6T*O 00 30061*0 00 30005*0- 00 50057*3 00 30057*0 

T 10-36SSa*D V0-3900f*0 60-50701 *0 TO 36i.8E*0 90-39611*0- 00 50006*0- 00 30005*0- TO 30001*0 

Q 

27 n XX kX n iu xx~ 

******************** SS3aiS »»*»»»*♦»♦»*♦»♦ ♦»♦»♦ *» »»♦♦»«»»» b3iN33 SS3diS ♦*•♦*****♦** 

i^1ttin,t*****lf*t¥**¥***********m** S3miM»ri3 SS331S 3MHnndD3 *¥»¥*¥*¥* ¥**¥*i'****ii*i***¥iT 


TO 36i.0E*O 10 30061*0 01 7 

TO 3610C*C TO BOOST'D 07 1 

SSTili dJ'n'53 *TM 

3Mi0 333 3 3AU.03333 J^3W313 


TO 35181 *0 TO 30671*0- 00 30679*0 00 3 0 679*0 00 306 16*0 00 30057*0- 00 50571 *D 00 30671*0 70 31511*0 10 31011*0 01 7 

90-35757*0- 10 30671 *0- TO 30671*0- 10 30057*0 10-3^86^*0 00 30057*0- 00 30067*0- 00 30006*0 70 3161X^0 — TO 31flt*0 07 T 


^ ^ ^ ^ 71 n iMlVlHitOD WlN3w'3dDNl *d*I 

— *••**•* **•*!* ••••* 3Al iyndlO »♦»»»*♦»»♦ »» p Itf lNad3 aPNI «««*#«*«0»*** ««** DliSVid ***** I'il3M3'l3 

♦♦♦♦•********************************** SNlYdlS OllSVId ************************************ “ '■■■ 


00 

30069*0 

00 

3000£*0- 

00 

3006€*0 

00 

3006£*0 

90-3891^*0 

90-37611*0 

90-36^£8*0- 

90-3>9£7*0 

TO 50507*0 

00 30055*0 

01 

7 

10 

-37187*0 

00 

3000e*D- 

00' 

300 0E*0‘- 

i6' 

'30001 *0 

90-'3Biii *0 ■ 

90-37611*0 

90-31C66*0 

90-30966*0 

TO' 30607*0 

00 aD05S*0 

07 

T 


AX 


12 


AA 


XX 

AX 

?Z 

AA 

XX 

3Ali?indflO 

“TflTN3rtIdO'Ni 

■■* J*T 

■■ •ON 


*»»»»»»♦**♦»»**♦*» BAIiVlOdaO ■»»■»»»»**»»» «. ♦»»»»»»»»»»»*T»** 1VlN3l>)3d JNjl ♦*»*♦»*»»* »** »♦* S NIVdiS iVrtdBHl *»* i'<3w313 

o********* ***************************** SNWaiS DI1SV13 **********************it*|i*********** ■ 



6 €-l ‘6 


INCFEMEM 


11 


CLw^LATIV E THEP ^AL LC ACS aNC Z-LCA CS FCS LCAC IKCPEKEKT 
13 


ELEMrM I.C. 20 10 

TEMP. 0.13000E C2 C.130C0E_02 

2-LCAC 0.20000F CC C.2C0CCF OC 


rTEPATIve ERRCP 

ITE PATIV E EPRCP 

fTFPATfVF FRPCP 
nfPAT IVE_.EPRCP 
ITEPATIVE EPRCR 
. .. . LTEPATI VE.iPPCP 
ITEPA7IV»^ EPRCP 

ITERAT IV E FPRCP 

ITERATIVE EPRCP 
_ ITEP^IVF_fRPCP 
ITEPATivE ERROR 

iterative errcp 

iTEPATIVF ERROR 


» 0.S238CE 00 

- O.EieOfeE 00 
= o.SfiewF 00 
= 0l5<555?E CO 

= 0.A3190E CO 

* . 0 ,iS 9 .gJlE::llL 

* O.32527E-01 
C.9CABeE“02 

= 0.27425E-02 

* 0.78IQ 2g-Q3 

= O.liTsiF-O'* 
= 0. 15091 E-05 

= 0.12363E-05 


t 


END CF LCAC I^CP£^E^T 13 


.INCREMEKT 13 

MECI-AMCAL LCAC CURVE FACTORS = C.75C0E Oli O.C ^ 

CRJF.P TIRE IKCPEMENT l._.OiC S 

NO. ELASTIC ELEMENTS ® C, NO. PLASTIC ELEMENTS =2 ' ® 

0 EL EMENTS HAVE CHANGEC ELASTIC TO PLASTIC. 0 ELEMENTS PLASTIC TG ELASTIC DLRING THIS INCREMENT 
SPECIFIFC MAX. NO. STIFFNESS UPDATES = 3. NC. UPDATES PERFCPMEC = 1 

SPECIFIFC RAX. NC._ IlfRATICNS PER UPCATE = 10« NO. ITERATI ONS PEPFQPME C SINCE LAST LPDATE = 3 

specified rax. liNBALANCEC-FCRCE EPPCP = O.IOCOE-04 t ACTUAL ERROR »' 0.'i236F-o'5 


*F**.*^* . C UMUL ATIVE INTERNAL FORCES AND DISPLACEMENTS ♦***» *♦* 


NODE 
Ls L 

** 

««*»**i**i. FORCES 
U V 

***** CISPLACFMENTS ****** 
U V 







1 

3C 

C.1788148E-06 

0.224S992F Cl 

“O.4540O76E-C7 

C.7500CCOE 01 







2 

40 

0.3175164E-06 

0.2249993G Cl 

0.200a019F CO 

C.75CCCC0E 01 







3 

. 10 

0, 10C0672E-C6 

-0.2249993F 01 

0.0 

0.0 







4 

20 

-0.5963975E-C6 

-0.2249993E Cl 

C.2CC00CCF CO 

0.0 









9.1-40 i 


*»*»♦»*<♦*♦♦»♦«♦«»*»♦»♦»*♦♦**♦♦»*♦♦*» ELASTIC S T (t4IN S »»»♦*» ♦»»*»»» »**» »»♦ »**♦<*»♦»* *»* ♦♦♦*** 

ElEME^T STRAINS *♦* **♦**♦♦♦♦**♦* !NCR^I^<^NTAL ♦♦••♦*♦*********♦ '*********«*** CUMtLATIVF ••***•♦*•*♦♦•♦♦*** 


NO. 

I.c. 

INCREMENTAL 

CUMULATIVE 

XX 

YY ll XY 

XX 


YV 

zz 

XY 


1 

20 

0.9500E 

CO 

0.2500E 

01 

-0.7153E-C6 

0.5364E-06 -0.7153E-06 -C.3990E-O7 

C.IOOOE 

01 

-0. 30C0E 

00 -C.300CE CO 

-0.U26E- 

07 

2 

10 

0.4500E 

CC 

0.25 COE 

01 

“C.6557E«C6 

-0.4172E-C6 -0.2980E-06 -0.59606-06 

0.3500E 

00 

0.3500E 

OC -C.30CCE CC 

0.65C0E 

00 


I 

I 


plastic strains 

eLFMCNT ♦**•* plastic ViCRK •♦♦♦ *•♦♦♦**4***** 1NCPE>«ENTAL ****»********4*44 4**4t**Ai*4**** CUMULATIVE ' ♦•♦444***4**44**4i* 


NC. 

I.C. 

incremental 

CUMULATIVE 

XX 


¥Y 

ZZ 


XY 

XX 


YV 

ZZ 

XV 

1 

20 

0.2094E 01 

0.1400E 02 

O.SOOOE 

00 

-0.2500E 

00 -0.2500E 

00 

0.399aE-07 

C.3000E 

Cl 

-0. 1500E 

Cl -C.15C0E 

01 -0.212SE-06 

2 

10 

0.2094F Cl 

0.1400E 02 

C.1250E 

CO 

0.125CE 

00 -C.2500E 

00 

C.3750E CO 

0.7500E 

CO 

0.75CCE 

00 -C.1500F 

Cl Q.2250E Cl 


***»»*<p»4*»444444»»44»4»*»»*4 CUMULA T IVE STRESS Q UANTIT I ES 

ELEMENT EFFECTIVE EFFECTIVE **4*4i**4***4 STRESS CENTER •#*♦•**♦*******♦ 4*4444**»****#» STRESS 444**4***4*4***4**** 

UC, CENTER STRESS U U LI O U H U ^ 

__ l_ 20 0.20_00f Cl Q.4 50 0E 01 0.1333E 01 -C.6667E 00 -0«6667E 00 -0,929 S E-07 C.4500E 01 -0,US^E-C5 C.5253E-C5 -0.3B99E-07 

Z lb 0.2000E 01 0.4500E 01 0,3333E CO C.3333E 00 -0.6667E 00 O.IOOOE 01 0,2250E Cl 0.2250E 01 -0.1525E-C5 0.22S0E 01 


ELEHEM ErP_ SUM I N C ? E^E NTA L TCTA L S URFACE ♦*** EFFECTIVE PLAS T IC STRAINS ♦»* ***** EFFECTIVE CREEP STRAINS *** * 

NO,' I,C. CODE CODE TEMPERATURE TEMPERATURE YIELD SIZE INCREMENTAL SUM INCR. CUMULATIVE INCREMENTAL SUM INC«. CUMULATIVE ' 

1 20 0 2 C.IOOOF 01 0.1300E 02 0.25COE 01 0.50CCE CO C.5000E 01 0.3QOQE 01 0.0 C.30C0E Cl O.IOCOE 01 

2 C.IOOOE 01 0.1300E 02 0.2500E 01 0.5000E 00 O.SOOOE 01 0,3000? 01 0.0 0.3000E Cl O.IOOOE 01 
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CIKULATIVE THER»'J»L LGACS «NC Z-LCAtS FOP LCAC tNCRE>'E^T 


IKCFE^'E^T 


CLEf'FM !.C, 2G 10 

TEMP. 0.14COOE C2 C.14 QOOE O: 

Z-LCAC O.C C.6 


ITERATIVE 
Ilf.f'ATIVE 
i TEPATI VE 
I TPPATIVF 
ITEPATIVE 
ITERATIVE 
ITERATIVE 

iterative 

ITERATIVE 
ITERATIVE 
I terati vl 
ITERATIVE 


ERROR » 
EP^CP *, 
ERROR '» 
EPPCP = 
ERROR = 
ERROR = 
ERROR = 
ERROR « 
ERROR = 
EPRCP 5 
ERROR '■= 
ERROR » 


0.7E£4SE OO 
0. C160BE 00 
O.RZRaiE OO 
Q.EE 2S7E CO 
0.25E04E 00 
0.3 9«A5E-^ 
'0*li2Al8E-02 
0.127826-02 
0.22749E-03 
0.4C927ErO:l 
0.ld6 59E*-05 
0.93568E-05 



END OF LCAC INCREMENT 14 


INCREMENT 14 

U, _HEC_hAK_ICAj^ LOAC CURVE FACT ORS « 0«8800E Pit 0»0 

CREEP TIME INCREMENT * 0.0 

i.__NC^EL_ASTrC EJJME^S_= Of NO. PIASTIC ELEMENTS = 2 

“■ d" ELEMENTS HAVE CHANGED ELASTIC TG PLASTIC, 0 ELEMENTS PLASTIC TC ELASTIC DURING THIS INCREMENT 

SPECIFIEC MAX. NO. STIFFNESS UPCATES - 3, NC. UPDATES PERFCRMEC = 1 

SPECIFIED MAX, NO, ITERATICNS PER UPCATE » 10, NC. ITEPATICNS PERFORMED SINCE LAST UPDATE * 2 

specified max. UNBALANCED-FCRCE EPPCP = C.lCCOE-04, ACTUAL ERROR = 0.9357E-05 



•* NODE ** 


c umulative INTERNAL FCPCES AND C I S FI ACE RENT S *♦*»*♦** 
*«*****«« FORCES »**♦* DISPLACEMENTS •***♦* 


1 30 -C.4112 72CE-Q5 _0 .^24991 LJ3 L 

2 4C Q.4693330E-05 0.2624992F 01 

3 _ 10 -0.169 89C6E -D 5 -0.2E24992E Cl 

4 2C C.iua296E“05 -0.2624992E 01 


0«5436ea7E- C6 C,87S99S9E Cl 

0. 1668930E-05 C.8799999E 01 

C.C OiC 

-0.6556511E-06 0.0 
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TD 30301*3 to 30D0C*3 O'O 10 3000^*0 lO 30309’0 73 30007*0 TO dDSi^'0 10 300^1*0 70 30007*0 e 0 07 2 

70 30007*0 to 3000E*0 0*0 70 lOOD^'O TO 30009*0 10 3000T*0 |{T30«l?*0 ZO 305^^7 ‘0 TO 30'001*0 2~~T5 3? T 

3Miyn^03 *ao\<i .^ns 7Vi^<3W3d3Ni d/Mitfionno -aoNi wns n?iN3»<3ii3Ni 37is aiaiA 3tiruyd3dw3i Sdfuyddddii 3Q30 aooo *o*i *on 
***♦ SNIVdlS d33d0 3 AU 33333 ***** *** SNiydiS OllSVTd 3A1133333 **♦* 3Dv3tfnS ItflJl lViN3M3d3NI WriS d-3 1N3W313 


CSJ 

I 

— liJ— - 

U 3 


CnJ 

70 

3S292'0 

90-3SieB’0 

- 70 35292*0 

10 

35292*0 

70 30527-0 

00 

3EEC8*0- 00 379Tt-*0 00 379Tb*0 

70 

30525*0 

70 

30052*3 

01 

2 

' i 
in 
tU‘ 

90 

-3tS92*0 

b0-38207“0 

S0-379bf0 

TO 

30525*0 

90-35907*0- 

00 

3Ce£8*0- 00 3Cet3'3- 10 37997*0 

70 

30525*0 

70 

30052*0 

02 

7 



AX 

77 

AA 


XX 

AX 


77 AA XX 


5S3dlS 


d3I>J30 

*3*1 

*3N 


*' 


I********** 

SSddiS ********«****«« 


3A1103333 

3AI 703333 



*•***•***«***«*««*««««*«*«**«**«« s3ij.iiNyrifi ss3dis 3AI itf irtf^rio ***************************** 


10 d^0€*0 10 30002*0- TO 3QOOT *0 tO 30007*0 00 30057*0 00 30005*0- 00 30052*0 OO 300SZ*0 20 37807*0 70 3579^*0 01 2 

90-3T5b«*0- i'O 30002*0- 70 30002*0- 7 0 3000V-0 90-39982*0- 00 30005*0- 00 30005*0- 70 30007*0 20 3780l*“O TO 3S78V*0 02‘ 7 


AX 77 n XX AX 72 AA 5(X 3ATl530dn0 7yi'J3W3dOMT *T*1 *W 

♦ ♦»»»♦ *»** **»*•»»* 3Miyinxr>0 *» ***♦♦» ***♦» **^********»***** 1ViN3>t3d3NI ************* **♦» >tdQH OliSVid ***** I'J3w313 


*************************************** SNlVdiS OliSVId ************************************ 



00 

30059*0 

00 

30DDE-0- OD 

3D0SE-0 

00 

3005C*0 

S0-3E1O7 *0- 

90-30965*0- 

90-30965*0 

90-37559*0 

70 

30082*0 

00 30ooe*o 

01 

2 


70 

-38999*0 

00 

300QE*0- DO 

300D£*0- 

10 

30007*0 

70-3^677*0 

90-30952*0 

90-30965*0 

0*0 

10 

30002*0 

03 3000€"0 

02 

T 



AX 


77 

AA 


XX 

AX 

77 

AA 

XX 

3AUyif7dn3 

iyiN3d3dON[ 

*3*1 

*0N 


«*******«*«**«***«__ 3AI jg3rt»<n0 *»»»*** *» ♦*** >**» »*»****♦♦*«** iyiN3*3d3Ml **»****«****♦ »** SNlVdlS 7yWd3H7 **♦ HBwaoS 

*•*******•***«****«*««««•*»****«*****«« SNIVdiS i)'llS9l3 ************************************ 


9.1-42 


e^-L*6 



ITERATIVE ERRCR » 0,^?388E-01 

ITERATIV E ERR OR « Q.<i6qsqE-Ql 

ITERATIVE e’rrcR « 0.2CSC6E-01 

QfJ-ATIVE ERROR « Q. IQ767E-01 

ITERATIVE ERROR * 0.A8A89E-02 

ITERATIVE E R RCP = 0.2<<238E-02 

ITERATIVE eVrCR = 0.37C40E-C3 

IJERATlVf ERR OR = Q.1B7L0£~03 

ITERATIVE ERROR * 0,7SC42E-04 

jurat IVE ER ROR = _ 0, 4 iaS4P-04 

iterative error » 0,f70COE-04 

ITERATIVE ERROR = 0.66C29E-0S 

ITERATIVE fPPCP ,= 0.SS181E>05 

_ _ _ ^ ^ _ 


JNCREREfia. _ i5 

MEf.I-AMCAL LCAC CURVE FACTORS = C.162SE 02t C.C 

CREEP TIRE INCREMENT = O.IOOOE 02 

NC. ELASTIC ELEMENTS = C» NG. PLASTIC ELEMENTS « 2 

Q_ ELEMENTS HAVE CMNGEO ELASTIC TO PLASTIC. 0 ELEMENTS PLASTIC TO ELASTI C DURIN G T HI S INCR EM ENT 

SPECIFIEC MAX. NC, STIFFNESS UPCATES = 3, NC. UPDATES PEPFCPMEC = 2 

SP ECI EIEC MAX. NO. I TER AT IONS PER UPDATE = IQ. NO. ITERATICNS PEPFCPMEC SINCE LA ST U PC4TE = 7 

SPECIFIEb MAX. UNeALANCEC-FCP'CE EPRCP * C.ICC0E-C4, ACTUAL ERROR = 0.S918E-0S 


44**_*** CU MULATIVE INTE R NAL F CPCE S AND DIS P LACE ME NTS «»»»» ♦♦». 
• * NODE 4* " *•*****■+*■ "fcPCES ********* ***** CISPLACEMENTS ****** 

NO. I.C. U V L V 

C.26822CSE-C5_ 0 129 99 967E C l 0. 84096 59E-06 0. 1E_2 874 9E 02. 

O.3903378B-C5 G.29'99991E 01 C. 1060135E-04 0.1628749E 02 

0,2552.fi89J-C5 ~0.2999989E CL C..0 _C,C 

0,3E54057F-b5 -3.2999990E 01 0.4021342E-05 0.0 


1 3C 

2 40 

3 1C 
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9.1-44 


************** *****»*****|»»^** * *** ** f LAST IC S TP/>IN S *^****»jM ♦♦^ *♦**•»•*•»♦*•♦>* ***♦♦**••* 

ELEM€^T ♦+>* THFPf^AL STRAINS *»♦ ************* INCBEKENTAL *****i»*********** 4****»*'**«<^i»« CLif^ULAT IVE *•*•'»•••»«♦•*»**•# 



. I. 

INCREMENTAL 

CUMULATIVE 

XX 


tv 

ll 


XV 


XX 


vv 

ll 


XV 


1 

20 

0.2363E Cl 

0.5162E 01 

C.IOOOE 

Cl 

-C.3000E 

00 -0.3000E 

00 

0.1A92E- 

C5 

C.20COE 

Cl 

-0.60COE 

CO -C,6CCCF 

CC 

C,1559C- 

C5 , 

2 

10 

0.2363E Cl 

0.5162E 01 

C.3500E 

00 

0.3500E 

00 -0.3000E 

00 

C*6500F 

GO 

C,70dOE 

CO 

0. 70CCE 

CC -C.60CCE 

CC 

o'.ncof 

Cl [ 


I 


»»**♦**»*♦♦♦»»♦*****♦♦** * *♦*»»***»»» PLASTIC ST R AIN S ♦»*♦♦*♦♦ ♦*♦*♦♦*»* * *♦» * »♦*»* ♦»»*♦♦»»»»»» 

EtEKENT ***** PLASTIC WORK *•♦* ' ****♦•*#***•* INCREMENTAL **•***♦*••**•*•** "*****»******« U L A T I V E *********Ts^*V**** 


NO. 


INCREMENTAL 

CUMULATIVE 

XX 


vv 

ll 


XV 

XX 


VV 

ZZ 


XV 

1 

2C 

0.5625E 01 

0.2450F 02 

O.IOOOE 

01 

-C.5000E 

00 “0.50CGE 

CO 

~C. 1347E-C5 

C.5000E 

Cl 

-0.25C0E 

01 -0.2500F 

01 

-0,18466-05 

2 

10 

0.562SE Oi 

0.2450E 02 

0.2500E 

00 

0.2500E 

00 -0.500CE 

CO 

0.750CE 00 

0.l25dE 

ci 

0-l2'50E 

ci -C.25CCE 

Cl 

0,37506 Cl 


*♦**♦*****»***♦♦»»»♦»♦♦»»♦♦♦»» *»*»»» _ CREEP strains ♦*•*****•***♦*•*••****••*******•••*•*** 

ELEMENT •***♦• CREEP WORK ***** ♦*•***♦**♦*♦* INCREMENTAL ****************» ***#*«**»*«** CUMULATIVE ****************** 


NO. 

I.C. 

INCREMENTAL 

CUMULATIVE 

XX 


vv 

ll 


XV 

XX 


VV 

ll 

XV 

LO_ _ _ 1. 

20 

0.I75BE 

C2 

0.2339E 02 

0.3125E 

Cl 

-C.1562E 

01 -0.1562E 

Cl 

-C.2273E-C5 

C.4125E 

01 

-0.2062E 

01 -0.20626 

01 -0.2994E-05 

1. 2 

10 

0.1758E 

C2 

0.233'9F 02 

C.7812E 

CO 

C.7812E 

00 -0.1562E 

01 

0.2344E Cl 

0.1031E 

01 

0.1C31E 

Cl -C-20626 

01 0,30946 01 













ELEMENT 

E-P 

SUM 

INCREMENTAL 

TOTAL 


SURFACE 

♦•** EFFECTIVE PLASTIC STRAINS **• 

**♦•• EFFECTIVE CREEP STRAINS ***• 

NC. I.C. 

1 

COCE 

0 

CCCF 
2. 

TEMPERATURE 
C-IOOOE 01 

TEMPERATURE 
Otl?COE C? 

VIELC SIZE 
0.3CCOE 01 

INCREMENTAL 
O.lCOOf Cl 

SUM INCR 
0. 7000E 

. cumulative 

01 0.5C00E 01 

INCREMENTAL SUM iNCR. CUMULATIVE 
C.3125E Cl C.6125E Cl Q.4125E CL 

2 10 

0 

2 

C.IOOOE 01 

0.1500E 

02 

0.30CCE Cl 

O.IOOOF Cl 

0, 7000E 

01 0.5000E 01 

0.3125E 01 0.ei25E CI 0.412SE Cl 



































•*^*♦*♦♦*♦*****♦*♦**^♦*****<1* CUMULATIVE stress QUANTITIES •**•♦****•************••*♦*•♦**•♦ 

ELEMENT effective" EFFECTIVE *»t********* STRESS CENTER **************** ' '******»i**i***^^. STRESS^ •»•**»*■*•*•■***♦*•*** 

NG, I.D. CfN TE P STPFSS )<X TV ll XY XX Y¥ ZZ XV 

1 20 0«3000E Cl 0.6000E 0 1 C.20QQE Cl -C.IOOOE CL - C«10C0 E C l -C.6373E~06 0.6000E, Cl -0. I130E-04 C.65A6E-05 0.355TE-C5 

2 10 0.3000E 01 0.6OCI0F 01 C.5CCCE CO C.500CE CO -0*1000f 01 C.15CCF 01 ""C.SOOOF Ci a^30C0F C 1 - C . 122CE- C4 C.30CCE Cl 
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9.2 THERMAL RATCHET 

This is a thermal ratchet problem, involving thermal cycling in conjunction 
with a sustained mechanical load. The finite-element idealization and 
mechanical loading are shown in Figure 9.2-1. The thermal loading consists 
of an alternate heating and cooling of the left half of the structure 
(elements 1 and 2). 

Because the stresses and thermal strains differ in the left and right 
halves of the structure, the BOPACE nodal -constraint option was used 
to allow vertical slip at the center. Thus the displacements at nodes 3-11 
and 4-12 are constrained to be equal in the X direction, but are allowed 
to have different values in the Y direction. Poisson's ratio is taken 
as 0.5 so as to avoid small errors which would otherwise be induced by 
intermediate yielding within an increment. 

Results are summarized in Table 9.2-1 for six increments, and the BOPACE 
input listing and printed output results are included at the end of 
this section. The mechanical loading is applied during the first 
increment and it then remains on the structure. In the second increment 
the thermal heating load is applied, and it results in plastic flow within 
the right side of the structure. Each succeeding heating and cooling cycle 
(two increments) results in continuing plastic flow and an increase of 
0.5 in displacement. Note that this occurs even though a part of the 
structure is always elastic, because yielding occurs during alternate 
increments in the left and right sides. This type of behavior must be 
considered in thermal -mechanical cycling of engines. 


9.2-1 
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PLATE = 2.0 X 1 .0 
THICKNESS = 10.0 
TENSILE YIELD POINT = 1 .0 
THERMAL COEFFICIENT 
OF EXPANSION = 1 .0 
E = 1 .0 
V = 0.5 


Y 



FIGURE 9.2-1: THERMAL RATCHET PROBLEM 


TABLE 9.2-1: THERMAL RATCHET DATA 


INCREMENT DISPLACEMENT TEMP. 1 & 2 TEMP. 11 & 12 


1 

0.75 

0 

2 

2,0 

1 .5 

3 

1.5 ^ 

0 

4 

2.5 

1.5 

5 

2.0 

0 

6 

3.0 

1.5 


9.2-2 


o o o o o o 






















SOLUTION HFT Hrn CODE = 3 

MAXIMUM NO, STIFrNFSS UPDATFS PFR INCREMENT = ? 

HA XI MUM total ITERATIONS PFR INCREMENT = LG 

MAXIMUM ELASTIC ITERATIONS PFR INCREMENT = 7 

MAXIMUM MAGNITUDE FOR ELASTIC-PLASTIC SUM CODE = 2 - 

MAXIMUM REDUCTIONS = I 

CON VERGENCE REOUCTIGN FA Q T OR = 0, S0003E 00 

MAx'i’mUM SPECIFIED ERROR NORM = 0,ioOOOE-OA- 

FRAC TION FROM END OF INCREMENT TO EVALUATE S L OPE = 0,IQOOOE 00 

PLANE-STRESS PROBLEM • 

NO, OF MATERIALS = 1 

OEFAULT THICKNESS = 0, LOOOQF 0 2 

FABRICATION TEMPERATURE = 0.0 




Ty TYP 


0 


MftTERIftL NO. 

1. temperature =0.0 

parade TEB. 

ISOTROPIC STRESS 

0*0 

<5, lOOOOE 03 

O.IOOOOE 01, 
O.IOOOOE 01 

PARAMETER 

KINEMATIC SHAPE 

0.0 

Otiooooe 03 

0.0 , 
0.0 

TEMPERATURE * 

0.0 - . . 

PARAMETER 

ISOTROPIC STRESS 

0.0 

O.IOOOOE 03 

0. lOOOOE 01 

O.IOOOOE 01 

PARAMETER 

KINEMATIC SHARE 

0.0 

0, 0 


O.IOOOOE 03 0,0 










9.2-7 
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_ _ _ CUMlit.flTTVF THERMAL LDAnS ANH pflR iMrpcy = SiT 

INCREMENT \ TCOLP) 


ELEMENT ! . C. 1 2 11 '? 



END DP L04n [^!C‘^FM"NT 1 







SM7dlS aJjcO 


*0 0 ‘G 

0*0 7TK 0 'd“ 

dAIiVTni^nO ‘•a'ji'd ->lt5 IV iN; d oJi'4 J 
SNIVdiS jll'jVla 3AIi3J; = 


j'^i ^ U.1il A adnigdddkgi dafUvb3<Jlm31 3000 30UD 'O*! ‘QN 


jjVadnS iViGi ^Vi^^dtN3a^NI■“Rns d^3 iNlWlTi 


40 -3?)ifao *0- 
4.0-34-764 *0 

0*0 

t0-374fe£*0 


0*0' 90-56aSI*0- 00 j0034*0- 
0 *0 ; 

u*G ‘;0-37jfc4‘0 00 30054*0 

0 '0 90-3ttU*0- 00 30034*0 


00 30054*0 


0*0 00 30064*0 

■0’’0 00 30054*0 


0*0 z 

i 5 *o‘^r 




ddii-jjQ 5 S 3 aj.d 

S 3 1 1 J iN vriO SS 3 bi S -i A 1 1 V ll'Wfi J 


3AI40d333 3A11033J3 IN3H313 





li aA XX AX 

r- aAli.? llU-i)3 ♦ y A * ;ji i,. ^ a 

yyjjOjsyyyy StlJlVblS jllSV 1 <J 


ay»yyyi(ty I 7 xN 3 W 3 b Jiv i y yi;!»y y 

vvy ...»! fiyyyyy *, yyyyyyy *.' vyy **** y # y » jjiy »*!* 


3AIi7lfmnD 1Vi.N3H3«3NI "O*! W 
**•# >bG(i JllStld ***♦♦ 1 N 3 M 3 "B 


90-30671*0- 00 3Gyic*C 


90-3 £ o 1 1 *0 
0*0 


00 jL)b4t.*C DD 30344*0 00 dCC54'0- OOTIt 

00 dOatc’C- 00 30b4t*0- 00 30054*0 


- 00 30544*0- DO 30054*0 


t *0- 00 30545*0 (.0 jObyt-'O CO d005i.*0- 













ELEMENT l.D. 1 2 11 12 

Q> 1 50_00E 01 0> 15000F 01 0. 0 Q^_Q 

z-Ln*o 0,0 0,0 0. 0 0,0 


ITERATIVE ERROR = 0, 12500E 01 

ITERATIVE ERROR = 0.62744E-05 

ITERATIVE ERROR = 0.3535SE 00 

HERA T IV E E RROR <= 0. 3535feE 0 0 

ITERATIVE ERROR = 0,35356E 00 

ITERATIVE ERROR := 0. ?337(1P 00 

ITERATIVE ERROR = O.Ul<e2E 00 
_ _ITERATIVE ERROR = 0.78574E>01 

ITERATIVE ERROR = 0,4l5<?9E-0l 

ITERATIVE ERROR = 0 . 21435E-01 

ITERATIVE ERROR = o’, 5'4369E-05 



Q, lOOQE 01, 0,0 


MECHANICAL LOAD CURVE FACTORS = 

CREEP TIME INCREMENT = 0,0 

NO, ELASTIC ELEMENTS = 2 , NO. PLASTIC ELEM ENTS = 

2 


ELEMENTS HAVE CHANGED ELASTIC TO PLASTIC, 
SPECIPIEO MAX, NO, STIFFNESS UPDATES = 2. NO. 
SPECIFIED MAX. 


0 EL FM ENT'S PLAS't'IC TO FL aSTTC DURING THIS iMrpry.PNjT 

PFR^TpRMFD = \ 

ND, ITERATIONS PEP UPDATE * 1 0 , ND, T ^ER ATION S p F= OP ‘-IF wr f= LAST UPDA TF = ? 

SPECIF lEO MAX. UNBAL ANCE D-FQRCE ERROR ~ 0. 10 ODE -04 . ACTUAL ERROR = 0.7'^?0r-n'^ 



** NODE ** *ai*3A***ii PqRCES **«ir***:*c± *3**.^ D I 5 P L A C F N >■ N T S 



U 


-0.4999997E 01 -0 .■'9e357Se-0S 


0. 0 


COOrtc, -'.c, q g;-| 



l-SQ 






ELEMENT *** THERMAL STRAINS 

INCREMENT AL CUMULATIVE 




ft * A # \\ 

X Y 


. I Q»1S00 E 0 1 0»1SQ0E 01 ^C»25QQg O'! 0. 1~>=^0" 00 00 

? 0.1500E 01 0.1500F 01 -Q,2500E 00 0. 1’SC= OQ 0. l^^'^OF 00 

1 1 QaO SiO --OtSSOOi Qii 0.1?=ior O.I??QE ^0 

12 0,0 0,0 -0.2500F QQ O.USOF 00 Ci.l?50'' 00 


: on - 0 . ^ “ ^ 00 ; o_o_* i._ 2 ?cci=' oo -o. 3is®E-06 

0.'!+76'PF-OA 
0 , q ? 3 Fc - o'T 

-o.iooor 01 o.^oooc no o.Fooor oo '' 0 . 202 '’f- 06 ’ 


******** ***#* Ki *** AA ^' 


element ***** PLASTIC WORK **** ****-********* ! 'C M *; 

1^LCR_£ME^IT AL_ CUMULATIVE XX YY f? 


0| 'FTIC STRAINS 


************* 


****************** 
7 XY 


0.0 


0.0 o.n 0.0 0,0 

-0. 1000*^ Q1 0. FOOiir 00 J. g^OOOE 00 0.7732^-07 

12 0,1000? 01 O.lOOOf 01 -0.1000? OV O.piO'^P? 00 n.FOOOE 00 - 0 . 15'=^OE-06 -O.IOOOE- 01 O.FOOOF 00 O.FOnOF no -0. 1F50F-06 


EFFECTIVE EFFECTIVE 
CFNTER_ ___ STRESS 


-1 I 0*0 

2 2 0,0 

^ U 

4 12 0,0 


****************** **** **** *** 
************ STRF SS C eWTc 
XX YY Z 


0, 50Q0E 00 Q.O 
0.5000E 00 0,0 
O.IOOOE 01 0.0 
O.IOOOE 01 n.Q 


C UM IJI. TIVF stress quant ITIFS t’^rc********.,******************** 


************ ,*PT * **#w*rt* ************ 

XX_ YY I 7 XY 


0 , SPOOF on -0.3I7 qf-o*> 0.0 

0.5000F 00 -0. 10?‘'E-0F. 0.0 

‘0, lOOOF 01 -O.F61feF-Ofe 0.0 
-O.IGOOF 01 0.?73SE-O7 , n.n' 


-0. 2126E-Q6 
0.3I7OF-06 

0.63S7F-07 

-0. i?piF-oe 


tl EHENT e-P SUM incremental TOTAL SURFiCF gFFECTlVF PLASTIC STRAINS 

NO, l.O, COOF CnOE TCMPERATUBE temperature YlFinSlZF T^!^.RFMF^'TA^ SUM UT=. '■ttmuiATrV'^ 


pperCTtvr CREEP STRAINS **** 
TNCFFMF^ITAl Si)m IKfrp, CUMUIATTVF 
0.0 0,0 















9.2-13 



CREEP TIhE increment = oTo ~~ ~ — ^ 

NQ . ELAST LC ELEMENTS ° 1 , NO. PLASTIC FI EMFNTS = 7 

2 ELEMENTS HAVE CHANGEO ElASTJC TT PLASTICt 2 Pt.ASTTC FLASTTC, ntjRtNG THIS TNCP'=M^^^T 

-SPECIFIE D MAX« N G.. STIFFNES'^ tlPDATFS = ?, NP. UPP^Tcc; pcpppoMpr) = 0 ' ' 

SPECIFIED MAX. ND» ITERATIONS PER UPDATE = lO» NO. I T FT A T in''! S pppcnt)H-n SINCE EAST UPOATF == c 

NrFD-EORCF FRRilR = 0. 1 OOPE-OA . AE Tj|Al FP.FOR = '). TOl a F-n^ 


^-^jL*** CUMULATIVE INTERNAL FORCES AND 0 1 c; dl AC EM ENT S 

NODE ** *i(t*****i** FORCES •***•■♦;*** ***** 0 I S PL AC E NTS 

— NQ« ItO. y y 
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ELEMENT *** THEffMaL STRAINS •*« 




**»****t**#*<(##**t^-ft***-^l(rfli, *5*. *■*«;*«, 


4t«*«**]^*«**«:* { Mr RFMFNTAL **na<-a~. 

YV 


FL A ST I ^ STRAINS ** *s<ii'*'****ite * ** * a ^<#** = * *** 


Ik It; 

XV 


C.500QP 00 - 0 »?FaOP oa_-0.?5aoc '•>(? 0.26?2I=-D6 Q.lOOOf 0! -Q. s nonf 00 -O.SOCG^^ on -0«SSTSP-07 

0.5000E 00 -0,?500E 00 -0. ^500? 00 -O.R<3A1E-07 O.IOOOF 01 -O.FDOOr 00 -0.5000E 00 0.3!17 aF-06 

S.5.Q00E 00 -O.?SO0£_OQ -a^2=^.Q0F 00 -0 . H ?0 g-06 -Q.SOOOr 00 0. R^OOF 00 0.7 500'^ 00 -0. i&iSiF-O? 


C.5000F 00 -C. 2F0DE 00 -0. ?50CE 00 0.?0'»6E-06 -0.5000F 00 0.7500F 00 0. 2^00= 00 0.59fiOF-Ofl 


ELEMENT ***** PLASTIC WORK •*«* 
0 . I , D ♦_ I NCR EMENT _Al_ CUM (JL A T I V E 


*******±*«*«±^**V!Att!k«^;\ 


I MC RFMFN 

XX YY 


I-***** A ST 1 C STRAINS *« *-**^ **>■»* 'c*>v *-*t s- -•;•!: t-i: * e a** *#*+^*-#*** 


I NC RFMF NT At **a****-,s;k'»*^***** ***** **-4 B CUMULATIV- ****************** 

YY 7 7 XY XX YV T» XY 


2 0.5000E 00 0,5000E 00 

JJ. .. Oj.0 Q.lOQOE 01 

12 0.0 O.IOOOE 01 


.SQOOE 00 -O.2S00F 00 -Q.R^OQP 00 -D . 1 2 7 0 . SPOOF 00 -0.2S00F 00 -3. ?SOOE 00 -0.5612E-07 

O.FOOOF 00 -0, 7S00F 00 -0. 2500F 00 0. r-07 O.SCOOT 00 -O. 7F00F 00 -0.7S00F 00 O.0941F-O7 

-0*0 OiO 0^0 0^0 -0. lOOOF 01 0. SnOQ!^ 00 O.'^OOCF 00 0.7732E-07 

0» 0 0.0 0.0 -O.iaOOF Ol O.fOCOF OO O.'^OOQF 00 -0. lSFOF-06 


EFFECTIVE 

EFFECT IVE 



0.0 

0*1 POPE 01 

0.0 

O.IOOOE 01 

Qai} 

M III! M ■ 1 

' 0,0 

0.5000E 00 


***************************** ru.NULATTVF STRFS 
************ STRFSS *********** «**** 


A NT I T t FS ****Jii »** ***!■: ***■<!***** ************ 



0, IDOO*" 01 0.317QF.-07 

0. 6000F 00 -0. ?Fft3F-07 


-0.5000F 00 O.1700F-O7 



0.2SP 3E-06 
0.1 109F-07 


0 
















INCREMENT 4 ( HOT ) 


0 Z-L[>.M.1S 


ELEMENT 1.0, 1 2 

11 

01 0. Q 0.0 

1? 







Z-LOAO 0.0 


0,0 , 

0.0 0,0 








ITERATIVE ERROR 

iterative error 

- 

0.125C0E 01 
0.36616E-05 









ITERATIVE ERROR 
. ITERATIVE ERROR 


0. 141A2E 00 
0.14143E 00 









ITERATIVE ERROR 
ITERATIVE ERROR 


0. 14143E 00 
0- loaTAE-05 









ITERATIVE ERROR 

= 

0.17333E-05 













E 

ND OF LOAD I 

IH 

■ 

Tl 

■- 





I.NC REMFNT <> < HOT ) 


MECHANICAL LOAD CURVE FACTORS = 
CREE P TIHE INCREMENT = Q.Q 


0. 1 oooe 01 f 0.0 


NO. ELASTIC ELEMENTS = 


2« NG. PLASTIC ELFMENTS = 


FLE^PNTS PLASTIC Tf.; ELASTIC OURINC THIS IMPPEMcnt 


~ SPECIFIED RAX, NO. STIFFNESS UPDATES » 

SPECIFIED MAX, NO. ITERATIONS PER UPOATF ^ 

2, NO. UPDATES o=P=DPMEP = 0 

10, NO. TPEPATJONS PER=nP«rO SINCE LAST UPDA = 

7 


fNJ 

SPECIFIED MAX, UNBALANCED-FORCE ERROR = C, 

lOCOF-04, ACTUAL EPFOP = O.t73?E-05 



(71 



♦**»♦** CUMULATIVE internal 

FORCES ANH 0 1 S ?L AC ^ M = NT S 




♦••NODE *• ♦♦*#***** forces **4<*4<>4*a 

NO. l.D. U U 

D I S P t. C FMENTS *•»».* 

U V - 





■ 0.0 0.0 





0.24999q9E 01 0 , 2029 1 76F-Q5 


O.ooo'-}9co= 00 
0. ?FOOonoF o: -o. ?o?oi?;r-os 


0. 2499997F 01 - 0. 1 PA77A 9E- 05 
Q. 499999^ 01 -0.1561245E-05 


0.2500001= 01 O.IOOOOOOF 01 
0.2500000= 01 -0. 29050! •' = -0'^ 


12 0.4999996F 01 0. 2 2960 ! 8 E“0 5 

13 -Q.499 0 997P 0 1 - 0. ? S709 Q6F- Q7 

lA' -0. 4999098E 01 -0 . 70656 ? 8 F-06 


O.250n0OL)E 01 
0 . 0 


0. 12500005 0! 
0. 0 


0,0 


o.i25on')0'‘ 0’ 


2-992£l-Sa 









eLEM6NT ♦** THERMAL STRAINS 
NO. _ l.D. incremental CUMULATIVE 


********1'** lk+**A* At A A A* tA aA 


A«A*****AAA*A INCRFMENTAL AASAtt 
XX YY 77 


AA * Att * AAA * AAA * A*A 

Y 


*****»[** a * A * 

XX 


* CtlMULST lYC 


A*** A A* A* A A*A 


A A n * AS * A * A *** AA**A 

: XY 


?SUO€ OQ 0 . 2 S 0 CF QO - 0 . •’T ^6 c-^;S O.SOnOE 00 - 0 , 2 ^ 00 ^ 00 - 0 . 2 ‘^OOE 00 - 0 . 3333 r -06 


2 0.1500E 01 0.1500E 01 -0.5000E 00 0.2SO0E 00 O.2SO0E 00 0.0 0. 5000F 00 -0.2FOOF 00 -0. ?S0C‘=' 00 0, 38 '^4 E-05 

11 O.Q OiO -0.50QQE 00 0.Z5OQE 00 0.250QE 00 -Q.^O^^Tr-Q? -O.lOOOr 01 O.^^OOOP 00 0. ^OOpr QQ -Q.4721F-07 

12 0,0 0.0 -0.5000E 00 O.2S0OE 00 0. ?500F 00 -0. 1792E-06 -0, lOOOF 0i O.SOOOF 00 0.^000'= 00 -0.1732E-06 


ELEHENT ♦***♦ PLASTIC WORK **** 
PEHENTAL 


****»***»*****»»*****»*^»»:***'^»<»-»'»J** PLASTIC strains At <• At a a A«: a Al'tC nt a AAA «Arta#1f! A A**# A**** 

♦•♦♦♦A******* incremental ****** aaa*a* *a# a* ************* CUMULATIVE 

XY XX YY 


Qi _0 Oi_g O.SOOOE 00 -C.g'^OOF 00 - 0 . 2800 E 00 • 0 . 5612 E -07 


.0 0,0 0.0 0.5000E 00 -0. ??00E OO -0. 2500E 00 0,R941F-07 

.?500E 00 0, 23_Q_0E 00 -O.lfrQRF-Q-^ -0,15Q0F 01 0,7S00E 00 0. 7SQCE 00 0.6173E-0? 
12 0,5000E 00 0,1500e 01 -0,S000F 00 0,2500E 00 0, 7500E 00 -O. ?043E-07 -0,1500E 01 0.7S00F 00 0.7S00E oO -0,1344E-06 


ELEMENT EFFECTIVE EFFECTIVE 
CENTER 


2 0.0 


12 0.0 


0.5000E 00 
0. lOQOF 


O.lOOOe 01 


»»**<t**»***»************a*A** cumulative stress QUANTITIES ******************************* *** 
************ STRESS CENTER **************** *************** $TRESS ********************* 

XV 




0 

0.5000E 00 

. Qa_Q 0* 0 Oa 


ksrmii 


0.50Q0E 00 -0,1‘’40F-05 0.0 

Q.IOOOE 01 >0.2Sg3E-0^ 0.0 


•O.lOOOf 01 0.343TF-06 ,0,0 


O.2503E-O6 

-0.3147E-07 


-0.11i5 5F-06 


SUH INCREMENTAL TOTAL 
I.O. CODE CODE TE.MPfRATURF TEMPERATURE 


ujvcaiiuatai 


SURFACE 

YIELD SIZE 




♦ *** cFPgr^TIVE PLASTIC STRAINS **=» 
,I NCRFMFNTAL sum incp, cumulative 


0. RrTOOF 00 0. 5000f 00 


***** effective creep strains **** 

INCPEMFNTAl SUM INCR. CUMULATIVE 

0.0 

















ELEMENT t.O. 

JEMP. 0.0 

Z-LOAO 0. 0 


I 2 

0.0 

0.0 


ITERATIVE ERRCR = 0.12500E 01 



INCREMENT 5 (COLOl 


■ i Wil ■■■±11^ ■ T M aifltlMa 


CREEP TIME increment = 0.0 


2 elements have changed ELASTIC TO PLASTIC, ? |=LEME^iT‘^ PLASTIC rn c-lASTIC DURING tmP', iNrPPMr^tT 

■SPECIFIED MAX. NO. STIFFNES S UPDATES = 2. NH. liPDATcs pc p c fic Mr r> = n 

SPECIFIFO MAX, NO, ITERATIONS PEP UPDATE = lO, N ^. TTFPATUIVS PFR'^HRmEO SINCE LAST UPDATf = 

ACtljal FPPna = O.apOSE-* 







T-;» ;■ !i;a *; ft ^ 

element *** THFRMAL strains *#* *•*«<********* ] ►!; oc Mf ji ( n^iv 

n 


cU'HJLiT JV = 


00 E 0 

1 0. 

0 


. 7RnOF 00 0. ■’n<:0r-06 O.IOOOF 01 -0,*- 


■ ■<**»* it y- 

XY 


0, 1000^^ 01 -0, ?00')r 00 -O. '^OOO'^ OO 0.3^’ftE-06 
q^SOTo!' 00 0 .?'^ 0 nr oO O.’tOOF. OO - 0 . 5136 F-OT 


- - V » ^ V - i v> ' ) ■ ) ui 'Tuvr. ; I ~w« 1 T or~ \j I 

0.5000E 00 -O.^'OO'^ 00 -0. 2500^ 00 0. ?0R6f^-0iS -0. 5000f 00 O.j'^OOF 00 0.2®^OOF 10 0,3'>2'»F-a7 


ELEMENT •>J>>i‘*t< PLASTIC WORK ♦*** 


NCRFMENTAL COHIH ATTVF 




«***-** I MC RpMR NT4(, 


***'^- ^ PL FC STRAINS 




* *»rt ■*** ************ 

XY 


O«lCiO0r 01 -O.FOCQP 00 - O.SOGOE 00 -0.77S3E-07 

0. lOOOF "01 -o.'^noor 00 -0. ‘^OOC'^ 00 0.1A53F-06 

-0»1F00C G! n.T^Qpr qq 0.7>=OOF OO 0. 617 3F-07 


- 0 , IfOO*^ 01 0.7 = ooe on 0 .?« 00 f= OO - 0 ,lRA 6 f- 0 f. 













IHCREMEHT 6 t HOT ) 


rUMULATfVE THERMAL LOAO*^ ANO Z-in^P<; IHAO 1^JCPF'■^F^IT f 


ELEMENT I,D. 1 2 

TEMP.. O.15QO0F 01 0.15000E 

11 

Qi 0-.0 0 . n 

1 7 



Z-LOAO 0, 0 


0.0 

0.0 0,0 




ITERATIVE ERROR 
ITERATIVE ERROR 

* 

0.12500E 01 
0. 39608E-05 





ITERATIVE ERROR 
ITERATIVE ERROR 

= 

0.14142E 00 
0.KU3E 00 





ITERATIVE ERROR 
ITERATIVE ERROR 

* 

0, 14143E 00 
3.12657E-05 





ITERATIVE ERROR 


0. IR410E-05 









E 

ND OF LOAD 1 

[NCR F H F T 

A 



_. IN CREMENT 6 t HOT I 

MECHANfCAL LOAO CURVE FACTORS = 0, 1 OOOE 01 1 0.0 

CREEP T I ME INCREMENT = Q.O ^ ^ 

NO, ELASTIC ELEMENTS = 2, NO, PLASTIC ELEMENTS = 2 ~ 

. 2 ELEMENTS HAVE CHANGE") ELASTIC Tn PLA STIC, 2 ELFM(-nTS PlASTlC ''H ELASTIC IXJ RI NG T « INCP.rMPNT 

SPECIFIED MAX, NO. STIFFNESS UPDATES = 2, NO, UPPATFS PER*=nPMFD = 0 

SPECIFIED MAX. NO, ITERATIONS PER UPDATF = 10, NO, ITFPATtnf^^ PEPt-nPMFD SINCE LAST UPDATE = 

M SPECIFIED MAX, UNflALANCED-FORCE ERROR = 0, lOOOE-OA, ACTUAL fRR'^P = 0. la-ilF-O? ' 
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AX i I AX 

SS-3dXS y y yif g i?.V iy s>i< v Jf>f i» sf tf »|i»» * * :f 

• J ( l «** lFl »# Ki «*; iA :,;*:? J * i ^ »».»=;! ** j *,| i ; V >,; i^JiliN^flR SS3^J.S 


U 




3AM33dd3 dAIlD3dd3 


_g0-J>lbI*0- U; aOOL'L’l) iO jpOOt'O TO j0Q0i^* 0- i.0 - 31 1 . £ T 'L« OU d0 0<»i '0 00 3 00‘s i; *0 00 3 000 S *0- 10 3 0002*0 00 30005*0 Z\ 

i0-35£i.£*0 ID 3000i*o [0 3L0qi*D ;0 aOoOZ'C- i.O-d it ^2’ c:- 00 30052*0 OC 30052*0 00 30 oCj 5 *3 fS 30003*0 U O 3O065 *'O 

90-jtiVl*0 00 30005*0- OQ 30 003*0- 10 30001 *0 0*0 0*0 0*0 0*0 TO 30001*0 

10-3tbi.i‘0- 00 3000i*0- 00 3 0 00 ,i * 0 - TO 30 001*0 0 *0 T’^U" 0”'^ 0^0 10 3000 1 *0 — ^ 


¥>? :iA1j.V mbfP ^ 

if Vir'i. ir^- 


iiifcrEmm 


=» vwag^ tsy: ■;;»»» i f :i.**’^* 3 V ti'J 3« d b OM «***♦*» >^*» ««44 OllSVTd 4**4* 1N3W3T9 

ivi#44^4 SNlVbib OliSVlrJ V=i- + »>p.4*444*:f.a*.4*»>.v4* *-» s» 191 44 4i» 4 v*4 4 4*4 


A0-3o2i.9*0- DO jOGOi'O 00 30005*0 TO 30 001 *0- VO-jItOl’C- 00 30052 *C 00 30052*0 

iO-3T*?£5*0- 00 30003*0 00 3D0oi*O TO 3000 T *0* '"iO-3 I *0- OU 300^2*0 00 300 52*0 

90-39 l£t.’0 00 30052 *0- 00 30Q5'J*0- 00 30 005 *0 0*0 00 dOOic'O 00 50052*0 

90-302t£*0- 00 5005d*0- 00 300ai‘O- GO 300 05*0 90- ao^-’O 2*0- 00 dODac *0 00 aOO 52 *0 


00 30052*0 00 30005*0- 

00 30005 *0- 
00 50052*0 00 30005*0- 


00 30052*0 DO 30005*0- 


0*0 

0*0 

10 3005 T*0 
10 30061*0 


0*0 21 

0*0 n £ 

TO 30051*0 Z Z 

10 30051*0 X T ■ 
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9.3 STAINLESS STEEL HYSTERESIS LOOPS 

This example is intended to aid the user in reducing cyclic test data into 
a form suitable for BOPACE input. The stainless steel specimen used in 
the example exhibits a pronounced variation in isotropic and kinematic 
hardening as the cycling progresses, and thus provides a severe test of 
the BOPACE capability for combined hardening. 

The cyclic test data are shown in Figure 9.3-1 by the solid lines.* The 
plotted points and dashed lines represent the analytical results obtained 
from a BOPACE run. The hysteresis loops are denoted by dash numbers, 
with the first number denoting the cycle number and the second denoting 
the 1st and 2nd half of the cycle. 

Because of the obvious variation in magnitude of kinematic hardening over 
the test cycles, the BOPACE option for variable kinematic hardening was 
employed in addition to the usual combined isotropic and kinematic hardening. 
The assumed hardening curves are given in Figure 9.3-2, and the BOPACE 
input data are listed at the end of this section. A summary of the analy- 
tical results is shown in Table 9.3-1. 

The Itest/analysis correlation is quite good, and could be further improved 
if one were willing to accept some amount of non-smoothness in the input 
hardening curves of Figure 9.3-2. An exact match would require hardening 
curves with discontinuous slopes at points corresponding to the strain 

* The test hysteresis loops were furnished by Dr. R, H, Mallett of the 
Advanced Reactor Division, Westinghouse Electric Corp. 


9.3-1 
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9.3 (Continued) 

range used in the cyclic test. This is not justifiable in general, 
because a proper test/analysis match would not result for other strain 
ranges. Although the test/analysis correlation is considered to be 
very satisfactory, cyclic results for other strain ranges would have 
to be compared before any definite conclusions can be drawn. 


9.3-2 
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TYPE 304 STAINLESS STEEL. ANNEALED 
TEMPERATURE = 1200°F 



STRAIN (%) 


FIGURE 9.3-1: CYCLIC TEST/ANAL YSIS CORRELATION FOR STAINLESS STEEL 


9.3-3 




KINEMATIC FACTOR, F KINEMATIC HARDENING, S (KSI) 





< » EFFECTIVE PLASTIC STRAIN (%) 
(C) KINEMATIC HARDENING FACTOR 


FIGURE 9*3-2: STAINLESS STEEL HARDENING ASSUMPTIONS 


9,3-4 




D5-1 7266-2 


TABLE 9.3-1: RESULTS FOR STAINLESS STEEL HARDENING 


3 


INCR . 

CYCLE 

DISPLACEMENT 

^XX 


1 1 

1 

0-2 

.05 

0.00 

8.90 

11.00 

2 

0-2 

.2 

0.70 

15.33 

14,64 

3 

0-2 

.4 

1 .74 

17.94 

• 16.20 

4 

0-2 

.7 

3.27 

20.44 

17.17 

5 

0-2 

1 .04 

4.97 

22.80 

17.82 

6 

1-1 

.9 

4.97 

- 2.12 

17.82 

7 

1-1 

.8 

3.38 

- 14.50 

17.88 

8 

1-1 

.6 

0.10 

- 18.02 

18.12 

9 

I-l 

.2 

- 2.66 

- 21.17 

18.50 

10 

1 -1 

-.4 

- 5.95 

- 25.01 

19.06 

11 

1-1 

- 1.04 

- 9.48 

- 28.97 

19.49 

12 

1-2 

-.9 

- 9.48 

- 4.05 

19.49 

13 

1-2 

-.8 

- 7.84 

11.65 

19.50 

14 

1-2 

-.6 

1 .05 

20.65 

19.60 

15 

1-2 

-.2 

5.89 

25.69 

• 19.82 

16 

1-2 

.4 

8.59 

28.70 

20.11 

17 

1-2 

1 .04 

10.83 

31,25 

20.42 

18 

2-1 

.9 

10.83 

6.33 

20.42 

19 

2-1 

.8 

9,85 

- 10.57 

20.42 

20 

2-1 

.6 

- 2,58 

- 23.07 

20.49 

21 

2-1 

.2 

- 8.64 

- 29.31 

20.67 

22 

2-1 

-.4 

- 10.54 

- 31.43 

20.89 

23 

• 2-1 

- 1 .04 

- 12.34 

- 33,45 

21.11 

24 

2-2 

-.9 

- 12.34 

- 8.53 

21 .11 

25 

2-2 

-.8 

- 12.06 

9.05 

21 .11 

26 

2-2 

-.6 

2.28 

23.43 

21 .15 

27 

2-2 

-.2 

9.10 

30.38 

21 .27 

28 

2-2 

.4 

11.00 

32.46 

21.46 

29 

2-2 

1 .04 

12.35 

33.97 

21.63 

30 

3-1 

.9 

12.35 

9.06 

- 21 .63 

31 

3-1 

.8 

12.35 

- 8.74 

21 .63 

32 

3-1 

.6 

- 2.80 

- 24.45 

21.66 

33 

3-1 

.2 

- 9.98 

- 31.73 

21 .76 

34 

3-1 

-.4 

- 11.71 

- 33.62 

21 .91 

35 

3-1 

- 1 .04 

- 13.09 

- 35.17 

22.08 

36 

3-2 

-.9 

- 13.09 

- 10.25 

22.08 

37 

3-2 

-.8 

- 13.09 

7.55 

22.08 

38 

3-2 

-.6 

2.25 

24.36 

22.11 

39 

3-2 

-.2 

9.83 

32,03 

22.20 

40 

3-2 

.4 

11 .49 

33.79 

22.30 

41 

3-2 

1 .04 

12.70 

35. TO 

22.40 


ITERATIONS 

1 

13 

8 

7 

7 

1 

9 

8 
7 
7 
7 

1 

7 
9 

8 
7 
7 

1 

7 
9 

8 
7 
7 

1 

5 

10 

7 

7 

7 

1 

1 

11 

8 
7 

7 

1 

1 

12 

8 
7 
7 


9.3-5 



M U» O 
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C/CLE 1-! 
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, CYCLF. 3-1. 

CYCLE ?-l 
. CYCLE 3-L , . 

"cycle 3-1 
CYCLE 3-1 
CYCLE 3- 1 T7 P 
...CYCL.E. 3-2. 

CYCLE 3-5 

C.YC.L E. 3-2 . 

CYCLE 3- 2 

_C_YU..t_3-2 

CYCLE 3^ rip 
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9,4 TYPICAL ENGINE ANALYSIS 

The BOPACE program has been directed toward thermo-structural analysis 
of thrust chamber liners in the Space Shuttle Main Engine (SSME) and 
in other scale model engines. This section describes a typical analysis 
for the SSME, 

Figure 9,4-1 shows a cross section of the SSME 470K thrust chamber. 

For the present demonstration a generalized plane-strain analysis is 
performed for a radial segment at a station 25.4 mm (1.0 inch) up- 
stream of the throat. Details of this segment are given in Figure 9.4-2. 
The finite element model is shown in Figure 9.4-3. Because this is a 
demonstration problem and because several loading cycles are analyzed, 
the selected mesh is fairly crude (only about 250 degrees of freedom). 
However the results appear to be quite good and indicate that a very 
fine mesh may not be needed for this type of analysis. 

The demonstration problem is summarized as follows. 

1. Three cycles at the given station in the SSME 470K engine are 
analyzed. 

2. Each cycle is defined by a start transient, sustained burn and 
shutdown. Sustained burn will last 500 seconds and shutdown is 
defined in the detailed report of the demonstration problem 
(Document D5-1 7266-3). 

3. For demonstration purposes, the following are included. 


9.4-1 
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9.4 (Continued) 

a. Creep behavior based on Boeing test data for NARLOY-Z. 

b. Estimated cyclic behavior of NARLOY-Z. 

c. Specified variation of z-strain from the inside of the liner 
to the outside of the structural jacket. 

d. Local tangential slip is allowed between the liner and Inconel 
jacket. 

A listing of the input data is included at the end of this section. In 
order to avoid unnecessary duplication, the listing includes thermal 
and z-load data for only the 1st increment (involving only thermal 
loads) and the 5th increment (involving both thermal and imposed z-strain 
loads ). 
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56 
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35 
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45 
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57 
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57 
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59 

55 

60 

92 
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56 

61 

93 
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61 

57 

62 

94 

1 

63 

59 

60 

95 

1 

64 
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61 
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1 

65 

61 

62 
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60 
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63 

98 

1 
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64 
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62 

66 
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64 

68 

63 

101 

1 

65 

69 
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66 
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67 

63 
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69 
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65 
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69 

70 
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117 
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80 

79 
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125 

1 

ftl 

85 

80 

126 

I 

82 

86 

31 

127 
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83 

79 
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2 

-.2286 

59 

2 

254 

63 

2 

-.254 

67 

2 

254 

71 

2 

254 

75 

2 

254 

79 

2 
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83 

2 

-.254 
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2 
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94 

1 

. 0762 

94 

2 

1524 

93 

1 

.1651 

92 

1 

. 1770 

91 

1 

. 0889 


20 


.31 1 . ,4 

.31 1 . .4 

.43 1,38 ,5 

I . 49 4.03 .45 

2 . 13 6 . 89 , 25 

2,98 9.65 .2 

4,15 13,44 ,2 

5 . 32 17 . 23 .2 

6,60 21.37 ,2 

8 . 31 26.89 .2 

10.01 32.41 .2 

II , 93 38.61 .2 

12 . 71 41.13 .2 

12.71 41.13 .2 
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12 . 71 41.13 ,2 
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12.71 41.13 .2 
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1 

257.65 

2 

257, 66 

3 

2 57 . 64 

L. 

257.63 

5 

257.61 

6 

257,62 

7 

257. 76 

8 

257.40 

9 

257, 40 

10 

257, 39 

11 

257.36 

12 

257.33 

13 

257.32 

14 

257.45 

15 

256 . 75 

1 6 

256. 73 

17 

256.68 

' 10 

256.64 

19 

256.60 

20 

256.57 

21 

256.64 

22 

256. 34 

23 

256 . 29 

24 

256.21 

25 

256.17 

26 

256.13 

27 

2 56.06 

20 

256. 13 

29 

255 . 43 

30 

255. 34 

31 

255.23 

32 

255.15 

33 

255.07 

34 

254.95 

35 

2 54 . 8 6 

36 

254.97 

37 

254. 80 

38 

254 . 75 

39 

254 . 6 -^ 

40 
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41 

254.41 

42 

254, 27 
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45 

253.67 

46 

253. 44 

47 

253. 30 

40 

253.08 

49 

252.96 

50 

253 . 30 
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253.09 

52 

252 .94 

53 
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54 

252 - 52 

55 

252. 32 

56 

252 . 13 

57 
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58 

251 . 78 
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251 . 51 
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62 
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2 50 . 51 
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65 

251. 30 

66 

250. 93 

67 

250. 46 

60 

250.00 
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71 

249 . 78 
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24 7 . 26 
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79 

247.33 
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83 

245 . 62 
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85 

244 . 16 
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244. 76 

87 
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92 
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239 . 97 
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2 37, 5 4 
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237. 75 
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107 

2 36.40 
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in 
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113 

2 34, 93 

114 
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U6 
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118 
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233. 76 

120 
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233.22 

122 

233.43 
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233.58 

124 

232.63 

125 

232. 92 
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23?. 97 

127 

232.30 
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232.55 
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232.69 
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231.92 
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232. 10 
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232.24 
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231. 92 
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2 32.2 4 
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2 31 ,44 
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231.60 
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2 31. 17 
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231.44 
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231,61 

142 

229.86 
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230.27 
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230.61 
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231. 27 
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231.47 
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230.42 
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231.08 
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231.30 
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230. 08 
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230.70 
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231. 12 
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231.39 
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230.38 
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230. S3 
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231.17 
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230.54 
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230.84 
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231. 09 
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23 0. 48 
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2 30. 76 

16 8 

230.96 

159 

230.65 
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230.90 
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2 30.9 5 
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230.75 
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230. 79 
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230.34 
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230.78 
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230.92 
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230. 82 

100 
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2 
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7 

311.78 
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3 09.27 

10 
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11 
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12 
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13 

308,92 

14 

309. 83 
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304.90 

17 
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18 

304.49 

19 

304,26 
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304.14 

21 
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302. 25 

23 
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301 .70 

25 

301.50 

26 

301. 28 

27 
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70 
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71 

252.39 

72 

255.05 

73 

25 3.26 

74 

250.74 

75 
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76 

251.79 

77 

246. 89 

78 

247.07 

79 

247. 70 

80 

246,91 

B1 

235.24 

82 

236. 39 

83 

2 37. 54 

84 
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85 
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86 
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87 

232.92 
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90 
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91 
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160. 29 

166 157,49 

167 1 50.78 

160 159.70 
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55 -0.000083 

56 -0, 000083 
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0. 000083 
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0,000085 
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0.000035 
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79 -0.000086 

80 -0.000085 



0,000087 
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C. 000089 
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0.000092 
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0.000102 
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Z-DIRECTION LOAD 
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0.000139 
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A,0 INPUTS INTERPOLATION AND DATA GENERATION 

Purpose - The program INPUTS was developed to serve as an interface 
between thermal analyzers such as BETA (Boeing Engineering Thermal 
Analyzer) and BOPACE. For the BETA program the structure is modeled 
as a network of lumped masses which in general do not correspond to the 
finite element nodes required for BOPACE. To expedite and automate 
preparation of thermal input data for BOPACE, INPUTB operates on the 

output from BETA to produce element temperatures for input to BOPACE. 

The generation of thermal data for BOPACE is shown schematically in 
Figure A. 0-1. Note that the BETA thermal analyzer can be replaced by 
any other thermal analyzer (e.g., SINDA, NASTRAN). 

However, INPUTB is of a general nature, and it can be used equally well 
for interpolation and data generation of the user-prescribed element z- 
di recti on loads. The INPUTB program will be discussed here in terms of 
temperature interpolation, but its use in interpolating other quantities 
should be obvious. 

Method - If temperatures are given as functions of time for some two- 
dimensional mesh of points, then the temperature at any other point and 
time can be determined by a spatial and temporal interpolation. 

The time variations of temperature are computed by a simple linear inter- 
polation. while the spatial variations are treated with an area weighting 
process equivalent to linear interpolation. 
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Figure A.0-1 : THERMAL ANALYSIS/GENERATION OF THERMAL DATA FOR BOPACE 
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FIGURE A. 0-2: SPATIAL INTERPOLATION 


DECK SETUP 



SELECTED BOPACE TIMES 
THERMAL NODE TEMPERATURE VS. TIME 

THERMAL NODE COORDINATES 
FINITE ELEMENT DEFINITIONS 

FINITE ELEMENT NODAL DATA 


FIGURE A. 0-3: INPUTS DECK SETUP 


A-3 


D5-17266>2 


A,0 (Continued) 

For the spatial interpolation, refer to Figure A. 0-2. Assume that the 
temperature is to be computed at point A. Temperatures at all other 
points are known. INPUTB logic proceeds as follows. 

(1) Find point 1, the point closest to A. 

(2) Locate point 2 as the next closest point such that the angle 
between lines A-1 and A-2 is greater than 90°. 

(3) Locate point 3 as the next closest point such that point A is 
wi thin the triangle 1-2-3. 

(4) The temperature at point A is determined from a weighting of 
the temperatures at points 1,2,3 according to areas A1 , A2, A3, 

i.e., 

(A1*T,M2*T2+A3*T3) 

(5) Special cases may arise where no triangle can be found according to 
the above procedure. These special cases have been accounted for. 

Summary of INPUTB Data - A pictorial of the INPUTB input deck is shown 
in Figure A. 0-3. The following is a listing of the input data by item 
(Formats are consistent with FORTRAN IV conventions). 

1. Finite Element Nodal Data (BOPACE Format) 

• For Each Node: Node I.D. 

I.D. of Coordinate System (0-Rectangul ar, 

1-Pol ar) 
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A.O (Continued) 


Nodal Coordinates (X»Y or R,e) 

I.D. of Displacement System (0-Rectangular, 
1-Polar) 

(215, 2F10. 0.15) 

• Blank Card after Last Node 

2. Finite Element Definitions (BOPACE Format) 

• For Each Element: Element I.D. 

Material Number 
Thickness 

3 Node Numbers (CCW Order) 

(215, FIO. 0,315) 

3. Thermal Node Coordinates 

• For Each Thermal Node: Node I.D., ICOORD, X(R), Y(6) 

ICOORD = 0 -Rectangular Coordinates 

ICOORD = 1 -R,0 Coordinates 
Up to 2 Nodes per Card 
2(215, 2F10.0) 

• Blank Card After Last Thermal Node 

4. Thermal Node Temperatures vs. Time 

• Number of Times, Default Temperature (I10,F10.0) 
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A.O 



(Continued) 

Time (FI 0.0) 

For Each Thermal Node: Node Number, Temperature 

Up to 4 Nodes Per Card 4(110, FIO.O) 

Blank Card after Last Node 

Repeat for each time up to number of times. 


For first time, default temperature will be assigned to 
unspecified nodes. For succeeding times, unspecified nodes will 
take on value at previous time. 


5. Selected Times 

For each time (BOPACE Time) at which Element Temperatures are to be 
Computed : 

• Time, ISTOP (FIO.O, 110) 

Note: ISTOP = 0, except on card following last time card 

set ISTOP = 9. 

Note: Max Input Times = 20 

Max Number of BOPACE Times = 12 

Max Number of Thermal Nodes = 300 (Max I.D. = 2000) 

Max Finite Element Nodes = 500 (Max I.D, = 2000) 

Max Number of Finite Elements = 800 (Max I.D. = 3000) 
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A.O (Continued) 

Output Data - Element Numbers and Temperatures are punched for each of 
the selected BOPACE times specified in the input as follows. 

• Time { F10.4 ,70(LH*) ) 

t Element Number, Element Temperature 4(110, FI 0.2) 
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